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Abstract. Using high resolution (R ~ 45 000), high S/N 
(~ 20-50) VLT/UVES data, we have analyzed the Lya 
forest of 3 QSOs in the neutral hydrogen (Hi) column 



density range iVm = 10 



12.5-16 



' 2 at 1.5 < z < 2.4. 



We combined our results with similar high-resolution, high 
S/N data in the literature at z > 2.4 to study the redshift 
evolution of the Lya forest at 1.5 < z < 4. We have ap- 
plied two types of analysis: the traditional Voigt profile 
fitting and statistics on the transmitted flux. The results 
from both analyses are in good agreement: 

1. The differential column density distribution function, 
f(Nni), of the Lya forest shows little evolution in 
the column density range A/hi = 10 12 - 5 ~ 14 cm~ 2 , 
f(N m ) oc JVjj-f , with /3 ~ 1.4-1.5 at 1.5 < z < 4 and 
with a possible increase of (3 to (3 ~ 1.7 at z < 1.8. A 
flattening of the power law slope at lower column den- 
sities at higher z can be attributed to more severe line 
blending. A deficiency of lines with TVri > 10 14 cm~ 2 
is more noticeable at lower z than at higher z. The one- 
point function and the two-point function of the flux 
confirm that strong lines do evolve faster than weak 
lines. 

2. The line number density per unit redshift, dn/dz, at 
A^hi = 10 13 ' 64-16 cm~ 2 is well fitted by a single power 
law, dn/dz oc (1 + z) 2 - 19±a - 27 , at 1.5 < 2 < 4. In com- 
bination with the HST results from the HST QSO ab- 
sorption line key project, the present data indicate that 
a flattening in the number density evolution occurs at 
z ~ 1.2. The line counts as a function of the filling fac- 
tor at the transmitted flux F in the range < F < 0.9 
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are constant in the interval 1.5 < z < 4. This suggests 
that the Hubble expansion is the main drive governing 
the forest evolution at z > 1.5 and that the meta- 
galactic UV background changes more slowly than a 
QSO-dominated background at z < 2. 

3. The observed cutoff Doppler parameter at the fixed 
column density A^hi = 10 13 5 cm~ 2 , 6 Cj i3.5, shows a 
weak increase with decreasing z, with a possible local 
b c ,i3.5 maximum at z ~ 2.9. 

4. The two-point velocity correlation function and the 
step optical depth correlation function show that the 
clustering strength increases as z decreases. 

5. The evolution of the mean Hi opacity, thi, is well ap- 
proximated by an empirical power law, thi oc (1 + 

z) 3.34±0.1 7j at 1.5 < 2 < 4. 

6. The baryon density, f2b, derived both from the mean 
H I opacity and from the one-point function of the flux 
is consistent with the hypothesis that most baryons 
(over 90%) reside in the forest at 1.5 < z < 4, with 
little change in the contribution to the density, O, as 
a function of z. 

Key words: Cosmology: observations - quasars: Lya 
forest - quasars: individual HE0515-4414, HE2217-2818, 
J2233-606, HS1946+7658, Q0302-003, Q0000-263 



1. Introduction 

The Lya forest imprinted in the spectra of high-z QSOs 
provides a unique and powerful tool to study the distribu- 
tion/evolution of baryonic matter and the physical status 
of the intergalactic medium (IGM) over a wide range of 
z up to z ~ 6. In addition, the Lya forest can be used to 
constrain cosmological parameters, such as the density pa- 
rameter Q and the baryon density f^, prov iding a test to 
current cosmological theories (Sargent et al. 198C ; Dave et 
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Table 1. Observation log 



QSO 



(mag) 



Wavelength 

(A) 



Exp. time 
(sec) 



Observing Date 



Comments 



HE0515-4414 

J2233-606 

J2233-606 

HE2217-2818 

HE2217-2818 



14.9 
17.5 



1.719 
2.238 



16.0 2.413 



3050-3860 
3050-3860 
3770-4980 
3050-3860 
3288-4522 



19000 
16200 
12300 
16200 
10800 



Dec. 14, 18, 1999 Reimers et al. ( |l99S| ) 
Oct. 8-12, 1999 
Oct. 10-16, 1999 

Oct. 5-6, 1999 Reimers et al. ( |l996| ) 
Sep. 27-28, 1999 



Taken from the SIMBAD astronomical database. 



al. |1999|; Im pey et al. |1999[ Schaye et al. [l999t Machacek 
et al. [2000| ). 

Although detections of C iv in the forest clouds suggest 
that the L ya forest is c losely related to galaxies (Cowie 
et al. 1995| ; Tytler et al. 1995| ), identification of its optical 
counterpart at z < 1 has produced different interpreta- 
tions: extended haloes of in tervening galaxies (Lanzetta 
et al. 1995 Chen et al. 1998 ) or Hi gas tracing the large- 
sca le dis tribution of galax ies and dark m atter (Morris et 
al. Il993t Shull et al. Il996l Le Brun et al. Il997|; Bowen et 



1998). Despite a lack of positive identification of opti- 



al. 

cal counterparts of the Lya forest, high resolution, high 
S/N data have provided a wealth of information on the 
cosmic evolution of the Lya forest, such as the space den- 
sity of absorbers, the distribution of column densities and 
Doppler widths, and the velocity correlation strengths (Lu 
et al. [I996t Cristiani et al. |1997| ; Kim et al. [1997D . 

Up to now, the systematic studies of the Lya forest at 
z < 1.7 have mostly relied on low resolution (R ~ 1 300- 
16 00 0) HS T observation s (Ba hcall et al. |1993| ; Weymann 
et al. 1998| ; Penton et al. 2000 ), which cannot be properly 
combined with the high-resolution (R ~ 40 000) ground- 
based data. Here, we present the observations of the Lya 
forest at 1.5 < z < 2.4 using the high resolution (R ~ 
45 000), high S/N (~ 20-50) VLT/UVES commissioning 
data on three QSOs. These observations take advan tage o f 
the high UV sensitivity of UVES (D'Odorico et al. |2000t ). 
Combining these data with similar Keck I/HIRES results 
at z > 2.5 from the literature, we address the z-evolution 
of the Lya forest at 1.5 < z < 4 as well as the physical 
properties of the Lya forest having Njji = io 125 ~ 16 cm~ 2 . 
In addition, when appropriate, we also compare our results 
with HST observations at z < 1.7. In Sect. 2, we describe 
the UVES observations and data reduction. In Sect. 3, we 
describe the conventional Voigt profile fitting technique 
and its application to the UVES spectra of the Lya forest 
in this study. In Sect. 4, we analyze the line sample ob- 
tained from the Voigt profile fitting. In Sect. 5, we show 
an analysis based on the transmitted flux or its optical 
depth, which supplements the Voigt profile fitting analy- 
sis in Sect. 4. This parallel analysis has the advantage of 
including absorptions with low optical depths which are 
usually excluded from the Voigt profile fitting analysis. It 



also gives a more robust comparison with numerical sim- 
ulations and other observations at similar S/N and reso- 
lution. We discuss our overall results in Sect. 6 and the 
conclusions are summarized in Sect. 7. In this study, all 
the quoted uncertainties are la errors. 

2. Observations and data reduction 

The data presented here were obtained during the Com- 
missioning I and II of UVES as a test of the instrument 
capability in the UV region and have been released by 
ESO for public use. The properties of the spectrograph 
and of its detectors are described in Dekker et al. (2000). 
Among the QSOs observed with UVES, we selected three 
QSOs for the Lya forest study at z ~ 2: HE0515-4414, 
J2233-606 and HE2217-2818. 

Complete wavelength coverage from the UV atmo- 
spheric cutoff ~ 3050 A to - 5000 A was obtained for 
the three QSOs with two setups which use dichroic beam 
splitters to feed the blue and red arm of the spectrograph 
in parallel. In this paper we discuss only the Lya forest 
observations which were recorded in the blue arm of the 
spectrograph. The pixel size in the direction of the disper- 
sion corresponds to 0.25 arcsec in the blue arm and the 
slit width was used typically 0.8-0.9 arcsec (the narrower 
slit being used for the brighter target HE0515-4414). The 
resolving power, as measured from several isolated Th-Ar 
lines distributed over the spectrum and extracted in the 
same way as the object spectra, is ~ 45 000 in the regions 
of interest. Table [I] lists the observation log and the mag- 
nitude of the observed QSOs. The exposure times are the 
sum of individual integrations ranging from 2 700 to 5 000 
seconds. 

The UVES data were reduced with the ESO- 
maintained MIDAS ECHELLE/UVES package. The indi- 
vidual frames were bias-subtracted and flat-fielded. The 
cosmic rays were flagged using a median filter. The 
sky-subtracted spectra were then optimally extracted, 
wavelength-calibrated, and merged. The wavelength cal- 
ibration was checked with the sky lines such as [O i] 
5577.3 38 A Nai 5989.953 A, and OH bands (Osterbrock 
1996 ). The typical uncertainty in wavelength is 



et al. 



O.OlX The wavelengths in the final spectra are vacuum 
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the complex. Although metal lines were searched for thor- 
oughly, it is possible that some unidentified metal lines are 
present in the Hi line lists from VPFIT. In most cases, 
absorption lines with b < 15 kms -1 can be attributed to 
metal lines (Rauch et al. 1997). In our line lists, these 



heliocentric. The individually reduced spectra were com- 
bined with weighting corresponding to their S/N and re- 
sampled with a 0.05 A bin. The S/N varies across the spec- 
trum, increasing towards longer wavelengths for a given 
instrumental configuration. The typical S/N per pixel is 
- 20-50 for HE0515-4414 at 3090-3260 A, ~ 25-40 for 
J2233-606 at 3400-3850 A and ~ 45-50 for HE2217-2818 
at 3550-4050 A. 

The combined spectra were then normalized locally 
using a 5th order polynomial fit. There is no optimal 
method to determine the real underlying continuum of 
high-z QSOs at wavelengths blueward of the Lya emission 
due to high numbers of Lya absorptions. The normaliza- 
tion of the spectra introduces the largest uncertainty in 
the study of weak forest lines. However, considering the 
high resolution of our data and the relatively low number 
density of the forest at z ~ 2, the continuum uncertainty 
should be considerably less than 10%. 

3. The Voigt profile fitting 

Conventionally, the Lya forest has been thought of as orig- 
inating in discrete clouds and has thus been analyzed as a 
collection of individual lines whose characteristics can be 
obtained by fitting the Voigt profiles. From the line fitting, 
three parameters are derived: the redshift of an absorption 
line, z, its Doppler parameter, b (if the line is broadened 
thermally, the b parameter gives the thermal temperature 
of a gas, 6 = y/2kT/m p , where T is the gas temperature, 
m p is the proton mass, and k is the Boltzmann constant), 
and its Hi column density, Am- 

We used Carswcll's VPFIT program (C arswell et al.: 
http:/ / www.ast.cam.ac.uk/~rfc/vpfit.htm] ) to fit the ab- 
sorption lines. For a selected wavelength region, VPFIT 
adjusts the initial guess solution to minimize the x 2 be- 
tween the data and the fit. We have chosen the reduced 
X 2 threshold for an acceptable fit to be 1.3 and we add 

more components if x 2 > 1.3. Even though the adopted forest (Songaila [l99|; Ellison et al. |1999[ ) 
threshold is somewhat arbitrary, a difference between the 
X 2 < 1.3 threshold and the x 2 < 1-1 threshold is negligi- 
ble when line blending is not severe, like at z ~ 2. Note 
that there is no unique solution for the Voigt profile fit- 



narrow lines are less than 5% of the absorption lines not 
identified as metal lines. Therefore, including these narrow 
lines does not change our conclusions significantly. 

We only consider the regions of a spectrum between 
the QSO's Lya and Ly/3 emission lines to avoid confusion 
with the Ly/3 forest. In addition, we exclude the regions 
close to the QSO's emission redshift to avoid the proximity 
effect. For HE0515-4414, we exclude a region of 4 500 km 
s" 1 , while for J2233-606 and HE2217-2818 we exclude a 
region of ~ 7000 km s _1 . 

Figs. 0, | and | show the spectra of HE2217-2818, 
J2233-606 and HE0515-4414, respectively, superposed 
with their Voigt profile fit (the fitted line lists [Tables A.l- 
A.3] from VPFIT with their errors are only electronically 
published) . The tick marks indicate the center of the lines 
fitted with VPFIT and the numbers above the bold tick 
marks indicate the number of the fitted line in the line 
lists, which starts from 0. From the residuals between the 
observed and the fitted spectra, the variation of S /N across 
the spectra is easily recognizable. Due to the limited S/N 
in the data, lines with Am < 10 12 ' 5 cm~ 2 become con- 
fused with noise. Therefore, we restricted our analysis to 
Ah i = 10 12 - 5 - 16 cm- 2 . We included all the forest lines 
regardless of the existence of associated metals because 
results from the HST observations are in general based on 
the equivalent width of the forest lines, not on the exis- 



tence of metal s (We ymann et al. 1998 ; Savaglio et al. 1999 ; 
Penton et al. 2000). In addition, the detection of metal 



lines in the forest at JVm > 10 cm and at z ~ 3 make 
it unclear whether there is a spread in the metallicity in 
the Lya forest or if there exists two different populations, 
such as a metal-free forest and a metal-contaminated Lya 



ting (cf. Kirkman & Tytler 1997). In particular, for high 
S /N data, absorption profiles show various degrees of de- 



parture from the Voigt profile (cf. Rauch 1996; Outram 
et al. 1999b| ). This departure can be fitted by adding one 
or two physically improbable narrow, weak lines, which 
results in overfitting of line profiles (see Sect. 5 for further 
discussion) . 

In this study, metal lines were excluded as follows: 
When isolated metal lines were identified, these portions of 
the spectrum were substituted by a mean normalized flux 
of 1 with noise similar to nearby spectral regions. When 
metal lines were embedded in a complex of Hi lines, the 
complex was fitted with Voigt profiles and the contribu- 
tion from the metal lines was subtracted from the profile of 



4. The Voigt profile analysis of the Lya forest 

In addition to the three QSOs observed with UVES, we 
have used the published line lists of three QSOs at higher 
redshift observed with similar resolution and S/N. Table || 
lists all the analyzed QSOs with their properties and the 
relevant references. We have avoided a region close to a 
damped Lya system in the spectrum of Q0000-263 and a 
lower S/N region in the spectrum of HS1946+7658. The 
spectrum of Q0302-003 does not include the region of the 
known void at z ~ 3.17 (cf. Dobrzycki & Bechtold 1991). 
The fitted line parameters with the associated errors of 
HS 1946+7658 and Q0000-263 were generated by VPFIT 
with a x 2 threshold of x 2 < 1-2 and x 2 < 1.1, respectively. 
The line list of Q0302-003 was generated by an automa- 
tized version of the Voigt profile fitting program by Hu et 
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Fig. 1. The spectrum of HE2217-2818 superposed with the Voigt profile fit. The residuals (the differences between 
the observed and the fitted flux) shown in the bottom part of each panel are shifted by —0.25. 
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Fig. 2. The spectrum of J2233-606 superposed with the Voigt profile fit. The residuals (the differences between the 
observed and the fitted flux) shown in the bottom part of each panel are shifted by —0.25. 
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HE0515-4414 
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Fig. 3. The spectrum of HE0515-4414 superposed with the Voigt profile fit. The residuals (the differences between 
the observed and the fitted flux) shown in the bottom part of each panel are shifted by —0.25. 



al. ( |l995|) with \ 2 ~ 1 an d the errors associated with the 
fitted parameters are not published. 

The results of profile fitting are known to be sensitive 
to the data quality as well as to the characteristics of the 
fitting program. As a consequence, comparing line lists 
obtained with different criteria is not usually straightfor- 
ward. Due to the use of a different fitting program, the 
line list of Q0302-003 at z ~ 2.9 should be treated with 
caution when combined with other line lists. A system- 
atic difference in b and AThi from VPFIT can introduce 
a slightly different behavior of the Lya forest at z ~ 2.9. 
While the difference would not change the study of the line 
number density or the correlation function significantly, it 
can affect the determination of a lower cutoff b envelope 
in the N^j-b diagrams. Furthermore, the six QSOs in Ta- 
ble a cover the Lyo; forest at 1.54 <c zj^y a < 4.0 with a 
fairly regular spacing. There is very little overlap between 
the Lya forests of the different QSOs and the effects of 
cosmic variance in the individual lines of sight might be 
important. 

4-1- The differential density distribution function 

The differential density distribution function, f(Nn t), is 
defined as the number of absorption lines per unit ab- 
sorption distance path and per unit column density as 



Table 2. Analyzed QSOs 



QSO 


-cm 


AA 






HE0515-4414 


1.719 


3090- 


-3260 


1.54- 


-1.68 


0.365 


J2233-606 b 


2.238 


3400- 


-3850 


1.80- 


-2.17 


1.104 


HE2217-2818 


2.413 


3550- 


-4050 


1.92- 


-2.33 


1.286 


HS1946+7658 C 


3.051 


4252- 


-4635 


2.50- 


-2.81 


1.157 


Q0302-003 d 


3.290 


4410- 


-5000 


2.63- 


-3.11 


1.878 


Q0000-263 6 


4.127 


5450- 


-6100 


3.48- 


-4.02 


2.540 



For go = 0. 

See also Cristiani fc D'O dorico ( |2000| ) 
Kirkman fc Tyt ler ( |L997| ) 
Hu et al. f|l995|) 
Lu et al. (jl996j) 



a function of A^hi (equivalent to the luminosity func- 
tion of galaxies). The absorption distance path X(z) is 
defined as X(z) = i[(l + z) 2 — 1] for q a — or as 
X(z) = + zfl 2 - 1] for q a = 0.5. We used q a = 
for dX to compare our /(A^hi) with the published /(ATm) 
from the literature (Table || lists the values of dX for 
Qo = 0). Empirically, /(ATm) is fitted by a power law: 



f(N Ul )=AN- 
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Fig. [| shows the observed log/(Am) as a function of 
log Am for different rcdshifts. The dotted line represents 
the incompleteness-corrected /(Am) at z ~ 2.8 from Hu 
et al. ( jl99i| , i.e. /(Am) = 4.9 x 10 7 7V Hl 146 . Note that 
the apparent flattening of the slope towards lower col- 
umn densities in the observed log Am-log /(Am) diagram 
is caused by line blending and limited S/N, i.e. incom- 
pleteness, which becomes more severe at higher z. For 
incom pleteness -corre cted /(Am ) at z > 2.5 (Hu et al. 
1995| ; Lu et al. |l996| Kim et al. |l997| ), this apparent flat- 



tening disappears. The incompleteness-corrected /(Am) 
at z ~ 2.8 over Am = 10 12 - 5 " 16 cm -2 is similar to the 
incom pleteness-corrected /(-/Vhi) at z ~ 3.7 (Lu et al. 
1996| ) and at z ~ 3.2 (Kim et al. 1997 ) over the same 



column density range. The amount of incompleteness ex- 
trapolated from at z > 2.5 (Hu et al. 1995 ; Lu et al. 1996 ; 
Kim et al. 1997 Kirkman & Tytler 1997 ) becomes negli- 
gible at z < 2.4 and we assume the observed /(/Vm) as 
representative of the actual /(Am) at z < 2.4. 

In the column density range Am = 10 12 - 5-14 cm~ 2 , 
the observed /(Am) at z < 2.4 is in good agree- 
ment for the different QSOs and also agrees with the 
incompleteness-corrected /(-/Vm) at 2.6 < z < 4.0. This 
suggests that there is very little evolution in /(Am) in 
the interval 1.5 < z < 4 for forest lines with Am = 
10 i2.5-i4 cm -2^ At Am = 10 14 ~ 16 cm " 2 , / (Am) shows 
differences at different z. Kim et al. (1997) noted that 
at lower z, /(Am) starts to deviate from a single power 
law for Am si 10 14 cm~ 2 and that the column density 
at which the deviation from a single power-law starts 
decreases as z decreases. The deviation from the single 
power-law in /(Am) is evident in Fig. |j. While the forest 
at z ~ 3.7 is still well approximated by a single power-law 



over Am = 10 



14-16 



the forest at z < 
> 10 141 



2.4 starts to 
cm~ 2 with a 



deviate from the power law at Am 
decreasing number of lines at Ah i 

Table || lists the parameters of a maximum-likelihood 
power-law fit to various column density ranges. These col- 
umn density ranges are selected for compariso n wit h the 
previous obse rvatio nal results of Kim et al. ( 1997 ) and 
Pcnton et al. ( 2000 ) and wi th sim ulations of Zhang et al. 
( |1998| ) and Machacek et al. ( [2000D . At z ~ 2.1, the slope /3 
is approximately 1.4 in the interval Am = 10 125 ~ 14 cm~ 2 
and 1.68 in the interval A Hl = 10 14 " 16 cm" 2 , i.e. the 
slope is steeper for higher column density clouds. At 
z ~ 1.61, the slopes j3 ~ 1.70-1.72 are steeper for both 
column density ranges. This indicates that the slope of 
/(Am) increases from z ~ 2.1 to z ~ 1.6. Assuming a 
curve of growth with b = 25 ± 5 km s _1 , Penton et al. 
( |2000D found that the slope of /(-/Vm) at z ~ 0.036 over 
Am = 10 12 - 5 - 14 cm" 2 and over A Hl = 10 14 ~ 16 cm" 2 
is j3 = 1.72 ± 0.06 and /3 = 1.43 ± 0.35, respectively. The 
slopes over Am = 10 12 5-14 cm~ 2 are steeper at z ~ 0.036 
and at z ~ 1.6 than at z > 1.8, and suggest that the 
incompleteness correction at z > 1.8 might be underes- 
timated or that the slope becomes intrinsically steeper 
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Fig. 4. The differential density distribution function as a 
function of log Ah i without the incompleteness correction 
due to line blending and limited S/N. The vertical error 
bars represent the lc Poisson errors. 

at z < 1.8. The slopes over A H: = 10 12 - 8 ~ 14 cm" 2 and 
Ah i = 10 14 " 16 cm -2 are in agreement with the ones found 
by Kim et al. ( |1997| ) at z ~ 2.3. However, our measure- 
ment of [3 = 1.48 ± 0.15 at z - 2.1 (only from HE2217- 
2818 and J2233-606) over Am = lO 13 - 1 " 14 - 5 cm" 2 is lower 
than the previous determination of f3 = 1.79 ± 0.10 over 
the same column density range at z ~ 1.85 by Kulkarni et 
al. ( jTggg ). 

While these observed f3 values at 1.5 < z < 2.4 can be 



obtained with semi-analytic models by Hui et al. (1997), 
they are lower than the values predicted from numerical 



simulations (Zhang et al. 1998; Machacek et al. 2000), by 
more than 2a. The slope depends on the amplitude of 
the power spectrum and models with less power produce 
steeper slopes. Thus, the steeper slopes from the simula- 



tions by Zhang et al. ( |1998[ ) and Machacek et al. ( |2000[ ) 
suggest that their index for the power spectrum, n p = 1, 
might be smaller than the actual index of the power spec- 
trum. 



4-2. The evolution of the line number density 

The line number density per unit redshift is defined as 
dn/dz = (dn/dz)o(l + z) 7 , where (dn/dz)o is the local 
comoving number density of the forest. For a non-evolving 
population in the standard Friedmann universe with the 
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Table 3. The power-law fit of the distribution functions, /(AThi) = AN } 



z 


TVhi = 10 12 - 5 - 


- 14 cm" 2 


7V HI = 10 14 -°- 


- 16 cm" 2 


JVhi = 10 12 ' 8 - 


- 14 3 cm" 2 


N m = 10 12 - 5 - 


- 16 cm~ 2 


logA 


P 


logA 





log A 




log A 





1.61 


10.65 ±0.86 


1.72 ±0.16 


10.38 ± 0.35 


1.70 ±0.42 


10.85 ±0.74 


1.74 ±0.21 


10.47 ±0.88 


1.71 ±0.10 


1.98 


6.19 ± 1.07 


1.35 ±0.09 


6.49 ± 0.59 


1.40 ± 0.21 


8.02 ± 1.00 


1.49 ±0.11 


8.66 ± 1.09 


1.54 ±0.05 


2.13 


6.66 ± 1.11 


1.38 ±0.08 


14.48 ± 0.65 


1.94 ±0.24 


7.80 ± 1.05 


1.46 ±0.10 


9.05 ± 1.14 


1.56 ± 0.05 


1.94 a 


6.97 ± 1.28 


1.41 ±0.06 


10.43 ±0.80 


1.68 ±0.15 


8.16 ± 1.21 


1.50 ±0.07 


9.02 ± 1.30 


1.57 ±0.03 



1.54 < z < 2.33. 



cosmological constant A = 0[], 7 = 1 and 0.5 for qo = and 
0.5, respectively. In practice, the measured 7 is dependent 
on the chosen column density interval, the redshift and 
the spectral resolution. Therefore, comparisons between 



individual studies are complicated (Kim et al. 1997). 
Fig. |^ shows the number density evolution of the Lya 



forest in the interval TVhi = 10 



cm . This range 



has been chosen to allow a comparison with the HST re- 
sults from the HST QSO absorption line Key Project at 
z < 1.5 from Weymann et al. (1998), for which a thresh- 
old in the equivalent width of 0.24 A was adopted. We 
assumed the conversion between the equivalent width and 
the column density to be iVm = 1.33 x l0 20 W/\lf, where 
W is the equivalent width in angstrom, Ao is the wave- 
length of Lya in angstrom, / is the oscillator strength of 
Lya (Cowie fc Song aila 1986). The value of the square 
(Penton et al. 2000 ) was estimated under the assump- 
tion of b — 25 km s _1 from the equivalent widths (corre- 
sponding to the column density range TVhi > 10 14 cm" 2 ) 
and is lower than the extrapolated dn/dz at z ~ 0.04 
from the Weymann et al. results, but within the error 
bar. Pentagons (Savaglio et al. 1999] from the line fit- 
ting analysis) also correspond to the column density range 
A^hi > 10 14 cm -2 . Note that Ahi from W depends on 
an assumed b parameter, resolution and S/N. Also note 
that including lines with ATm > 10 16 cm" 2 from line fit- 
ting analyses introduces a further uncertainty on the line 
counting since different programs deblend completely sat- 
urated lines differently, resulting in different numbers of 
lines for the same saturated lines. 

The long-dashed line is the maximum-likelihood fit to 
the UVES and the HIRES data at z > 1.5: dn/dz = 
(9.06 ± 0.40) (1 + 2)2-19+0.27^ Thig ^ ig lo wer th an pre _ 

viously reported 7 ~ 2.75 ± 0.30 (Lu et al. |l99l[ Kim et 



1 Recent measurements from high-z supernov ae fa vor the 
non-zero cosmological constant (Perlmutter et al. 1999). When 



we use the mass density f2 m ~ 0.3 and the cosmological con- 
stant energy density £l\ ~ 0.7 for the flat universe as the re - 
sults from the supernova study favor (Perlmutter et al. 1999] ), 
the line number density for the non-evolving forest can still 
be approximated by a single power-law dn/dz oc (1 + z)~ ' 71 
at < z < 4.5 with a slight steepening at z < 1: At z > 1, 
dn/dz oc (1 ± z)~ ' 59 , while z < 1, dn/dz oc (1 + z)^ 1 ' 15 . 



al. 1997). This slope is steeper than the expected values 
for the non-evolving forest for a universe with £l\ — 0.7, 
£l m = 0.3 and £1 = 1. These results suggest that the Lya 
forest at TVhi = 10 1364 - 16 cm- 2 evolves and that its evo- 
lution slows down as z decreases. Interestingly, the HST 
data point at < z > = 1.6 (the open triangle at the bound- 
ary of the shaded area) , which has been measured in the 
line-of-sight to the QSO UM 18 an d was suggested to be 
an outlier by Weymann et al. ( 1998 ), is now in good agree- 
ment with the extrapolated fit from higher z. 

Despite the different line counting methods between 
the HST observations (based on the equivalent width) 
and the high-resolution observations (based on the profile 
fitting), a change of the slope in the Lya number den- 
sity does seem to be real. The UVES observations suggest 
that the slow-down in the evolution does occur at z ~ 1.2, 
rather than at z ~ 1.7 as previously suggested (Impey 



et al. 1996; Weymann et al. 1998), although the different 
methods of line counting at higher and lower z make it 
a little uncertain. At least, down to z ~ 1.5, the number 
density of the forest evolves as at higher z, which sug- 
gests that any major drive governing the forest evolution 
at z > 2 continues to dominate the forest evolution down 
to z ~ 1.5. Since the Hubble expansion is the main drive 
gover ning t he forest evolution at z > 2 (Miralda-Escude 
et al. 1996 ), the continuously decreasing number density 
of the forest down to z ~ 1.5 implies that the Hubble ex- 
pansion continues to dominate the forest evolution down 
to z ~ 1.5. 



Fig. H is similar to Fig. ||, except for the A^hi range: 
A^hi = 10 131 - 14 cm- 2 . The correction for incomplete- 
ness due to line blending is still negligible in this col- 
umn density range (Fig. |] shows that the number of lines 
per unit column density over iVm = 10 131 - 14 cm -2 is 
still well represented by a single power-law). Again, the 
square from Penton et al. ( 2000] ) is estimated from the 
equivalent widths with the assumed b = 25 km s . 
The dot-dashed line is the maximum-likelihood fit to 
the lower column density forest of the UVES and the 
HIRES data: dn/dz = (55.91 ± 2.00) (1 ± z ) 110 ±0-2i_ At 
2.4 < z < 4 and at 2.1 < z < 4, 7 = 0.90 ± 0.29 and 
7 = 1.00 ± 0.22, respectively. For the column density 
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Fig. 5. The number density evolution of the Lya forest. 
The column density range iVm = 10 13 - 64 " 16 cm" 2 has 
been chosen to allow a com parison with the HST results 
from Weymann ct al. ( 1998 ), which are shown as open tri- 
angles. Filled symbols are estimated from HE0515-4414 
at < z > = 1.61, from J2233-606 at < z > = 1.98 and 
from HE2217-2818 at < z > = 2.13, respectively. The 
star, open circles, and the diamond are taken from the 
HIRES dat a at similar resolutions by Lu et al . ( 1996|) , 
Kim et al. ( |l997| ), and Kirkman & Tytler ( |l997| ), respec- 



tively. The square and pentagons are taken from Penton 
et al. ( [2000D and Savaglio et al. (|l999| ), respectively, over 
AThi > 10 14 cm" 2 . Horizontal solid lines represent the 
z interval over which the number density was estimated. 
Vertical solid lines represent the la Poisson errors. The 
shaded area is the z range where UVES is extremely sen- 
sitive. The long-dashed line is the maximum likelihood fit 
to the UVES and the HIRES data at z > 1.5. The UVES 
observations indicate that the number density evolution 
of the Lya forest at z > 2.4 continues at least down to 
z ~ 1.5 and that a slope change occurs at z ~ 1.2. 



Fig. 6. The number density evolution of the Lya forest at 
the Njh range AT H i = 10 131 " 140 cm" 2 . Symbols are the 
same as in Fig. ||. The dot-dashed line (solid line) is the 
maximum likelihood fit to the lower column density forest 
including the < z > = 2.66 forest (excluding the < z > = 
2.66 forest), while the dashed line represents the maximum 
likelihood fit to the forest at AThi = 10 
Fig. | 



13.64-16 



from 



fitting, which, as discussed in Sect. 3, may occur especially 
in high S/N data. 

As previously noticed (Kim et al. 1997), the lower 



column density forest evolves at a slower rate than the 
higher column density forest. The evolutionary rate oc 
(1 + z ) 110±0 - 21 would be consistent with no evolution for 
q Q = 1 or mild evolution for q Q = 0.5. For = 0.7, 
f2 m = 0.3 and fl = 1, the Lya forest with A^hi = 
10 131 " 14 cm" 2 is mildly evolving at z > 1.5. The Lya 
forest with A^hi = 10 131 " 14 cm" 2 appears more numer- 
ous at z ~ than expected when extrapolating from the 
z > 1.5 range. 



range A^hi = 10 131 1 cm 2 , the forest does not show 
any strong evolution. 

Note that the point at < z > — 2.66 (diamond) from 
Kirkman & Tytler (1997) indicates a number density 
twice as large as than at < z > = 2.87 in the interval 
jV H] = io 131 - 14 ° cm" 2 (excluding the < z > = 2.66 for- 
est, the maximum-likelihood fit becomes dn/dz = (47.77± 
1.84) (l + z) 118±0 - 22 ). Although this discrepancy could re- 
sult from a real cosmic variance of the number density 
from sightlinc to sightline, the number density in the in- 



terval A^hi = 10 



13.64-16.0 



cm from the same line of sight 



is in good agreement with other HIRES data. The differ- 
ential density distribution function (Fig. ^) and the mean 
H I opacity (Fig. |l5|) towards this line of sight suggest that 
the discrepancy at = 10 131_ 



4-3. The lower cutoff b in the Nui-b distribution 
4.3.1. The lower cutoff b parameter 

For a photoionized gas, a temperature-density relation ex- 
ists, i.e. the equation of state: T — T (1 + §) 1T ~ l , where 
T is the gas temperature, To is the gas temperature at the 
mean gas density, 8 is the baryon overdensity, (pt — ~Pb)l~Pb 
(~p b is the mean baryon density), and jt is a constant 
which depends on the reionization history (Hui & Gncdin 
|1997 ). For an abrupt reionization at z > 5, the tempera- 
ture of the mean gas density decreases as z decreases af- 
ter the reionization, eventually approaching an asymptotic 
To. For a generally assumed QSO-dominated UV back- 
ground with a sudden turn-on of QSOs at 5 < z < 10, 
T decreases as z decreases at 2 < z < 4 (Hui & Gnedin 



19971) 
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Under the assumption that there are some lines which 
are broadened primarily by the thermal motion at any 
given column density, this equation of state translates into 
a lower cutoff 6(Am) envelope in the Am-& distribution: 
T and 5 can be derived from b and Am- For the equation 
of state T = T (1 + <5)^"\ 6 c (A H i) becomes 



log(6 c ) = log(6o) + (It - 1) log(A Hl ), 



(1) 



where log(6o) is the intercept of the cutoff in the log(-ZVni)- 
log b diagram and (IV — 1) is the slope of the cutoff (Schaye 
et al. 1999| ). The cutoff slope (Tt — 1) is proportional to 
(7t — !)■ This cutoff envelope provides a probe of the 
gas temperature of the IGM at a given z, thus giving a 
powerful constraint on the thermal history of the IGM (Hu 
et al 



1995; Lu et al. 199C; Kim et al. 1997 



Tytler p_997t Zhang et al. [1997| ; Scha ye et a l 

al " 



Kirkman & 
1999| Bryan 



2000; Ricotti et al. 



fc M achacek |2000| ; McD onald et 
2000| ; Schaye et al. |200C| ). 

In practice, defining 6 c (Ahi) in an objective manner is 
not trivial due to the finite number of available absorption 
lines, sightline-to-sightline cosmic variances, limited S/N, 
and unidentified metal lines. Among several methods pro- 
posed to derive & c (Ahi), we have adopted the following- 
three: the iterative power-law fit, the power-law fit to the 
smoothed b distribution, and the b distribution. We re- 
fer the reader to other papers for more methods to derive 
6 c (A Hl ) (Hu e t al. |1995| ; McDonald et al. |2000t Theuns & 



Zaroubi |200C| ) 



In our analysis, we divide the data points into 2 groups: 
Sample A and Sample B. Sample A consists of the lines 
in the range Am = lO 125-145 cm -2 with errors less than 
25% in both Am and b in order to avoid ill- fitted val- 
ues from VPFIT. Sample B consists of all the lines with 
Am = 10 12 - 5-14 - 5 cm" 2 regardless of errors. The criteria 
for Sample A and Sample B are chosen to compare our 
results with the previous results by Schaye et al. (2000) 



and to investigate whether it is reasonable to include the 
Q0302-003 line list for which error estimates are not given. 
As no errors are available for Q0302-003, no Sample A 
can be defined at < z > = 2.87. Note that including the 



relatively few lines with Am = 10 



14.5-16 



cm 2 does not 



change the results significantly. There is hardly any over- 
lap in z, except for J2233-606 and HE2217-2818. Since one 
of our aims is to probe the z-evolution of 6 c (Am), we an- 
alyze the b distribution of each line of sight individually 
to derive 6 c (Am). 

4.3.2. The iterative power law fit 

Since the equation of state is a power law, it is reason- 
able to fit Ahi-6 c with one. We did so, using the boot- 
strap method described by Schaye et al. ( 1999| ), iterat- 
ing until convergence was reached. After each iteration, 
those points were excluded that lay more than one mean 
absolute deviation above the fit. Finally, the lines more 
than one mean absolute deviation below the fit were also 



taken out and the final power law fit, b c = cq iP A^ -1 , was 
carried out. The procedure was repeated over 200 boot- 
strap realizations in order to determine the full probability 
distribution of the parameters of the cutoff. As noted by 
Schaye et al. ( 2000D , the power law fit requires over 200 
available lines to reach stable fit parameters. 

Fig. ^ shows the iterative power law fit in the log Ah i~ 
b distributions. The noticeable difference between Sample 
A (cross symbols) and Sample B (cross symbols and open 
circles) occurs at Am < I0 13 cm" 2 . These lines usually 
come from blends or from weak, asymmetric absorption 
lines. Table || lists the fitted parameters, such as co, p and 
(Tt — 1), including b c values at the fixed column density 
Am = fO 13 5 cm" 2 , 6 c ,i3.5, for Sample A and Sample B. 

The power law fit between Sample A and Sample B 
docs not give a significant difference except at <z> = 2.13 
and at < z > — 3.75, for which several lines with b < 20 
km s" 1 and Am < I0 13 cm contribute to a different 
power law fit for Sample B. This suggests that using the 
Q0302-003 line list at < z > = 2.87 without error bars 
does not severely distort our conclusions. Note that the 
power law fit at < z > = 2.87 might be less steep with a 
higher intercept, if the same general behavior of errors also 
occurs for the Q0302-003 forest (larger errors at b < 20 
km s" 1 or 6 > 40 km s" 1 ). Due to the small number 
of lines (47 lines for Sample A and 56 lines for Sample 
B) at < z > = 1.61, the power law fit should be taken 
as an upper limit on & c (Am) and indeed it provides the 
highest 6 c .i3.5 among all the z bins. For both Sample A 
and Sample B, there is a weak trend of increasing 6 Ci i3.5 as 
z decreases, except at < z > = 2.87 which shows a higher 
&c,i3.5 value than at the adjacent z ranges (see Sect. 6.2 
for further discussion). On the other hand, the power law 
slope (Ft — 1) is rather ill-defined with z with a possible 
flatter slope at <z> = 3.75 than at z < 3.1. 

Fig. H shows the power law fit to Sample A at z ~ 
2.1 (small filled circles; 242 lines from J2233-606 and 
HE2217-28118) and at < z > = 3.75 (open squares; 209 
lines) over Am = f0 



12.5-14.5 



cm (upper panel) and over 
Am = 10 13 " 145 cm" 2 (lower panel). The fitted parame- 
ters are given in Table [|. For both Am ranges, the slopes 
of & c (Ahi) are steeper at z ~ 2.1 than at < z > = 3.75. 
This result, however, is certainly biased by the lack of lines 
with b < 15 km s" 1 and Am < 10 13 4 cm" 2 at higher z, 
due to the severe line blending. 



4.3.3. The power law fit to the smoothed b distribution 



Bryan & Machacek ( |2000D presented a method to measure 
6 c (JVhi) from a power law fit to a smoothed b distribu- 
tion, sorting absorption lines by JVm and then dividing 
them into groups containing similar numbers of lines. The 
b distribution in each group was then smoothed with a 
Gaussian filter with a smoothing constant 



exp(-(6i - b) 2 /2a 2 b ), 



(2) 
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Fig. 7. The logTVui-fr diagrams with the power law fits over iVm = l() 12 - 5_14 - 5 cm -2 . Errors are not displayed. Crosses 
indicate the data points with errors less than 25% both in iVm and b (Sample A), while open circles indicate the lines 
with errors greater than 25 % (Sample B consists of crosses and open circles together). At <z> — 2.87, no Sample A 
can be defined due to the lack of errors in the line list. Thus, Sample B was used as a substitute for Sample A when 
compared to Sample A from QSOs at other z. The solid line, the dot-dashed line and the dashed line at each panel 
indicate the power law fit to Sample A, to Sample B and to the <z> = 3.75 forest for comparisons, respectively. At 
< z > — 2.87 the solid line is not present since no Sample A can be defined. 



where Sbj (b) is the smoothed density of lines in each group 
j and i indicates the lines in the group. Then, the location 
of the first peak in the derivative of Sb,j(b) defines the 
lower cutoff at the average column density, Nui,j, for the 
j-th group. 

Fig. ^ shows the logiVni-6 diagram at each z with 
the b c (Nni) points for each group (filled circles) measured 
from the smoothed b distributions. We use the smoothing 
constant — 3 km s _1 . However, Sb,j(b) is largely in- 
sensitive to the smoothing constant. In general, 30 lines 
were included in each group except for the last group at 
higher iVm for which typically smaller numbers of lines 
were available. For this same reason, at < z > = 1.61 



groups of ~ 16 lines were used. The solid line represents 
the robust least-squares power law fit to filled circles: 
b c (Nni) = Co )S iV H I j r_1 . Table || lists the parameters of 
the power law fit to the smoothed b distributions. 

We find that the power law fit to the smoothed b distri- 
bution produces in general a lower intercept and a steeper 
slope than the iterative power law fit. It also produces 
smaller 6 c .i3.5 values. Direct comparison of Fig. |^ with 
Fig. [?] indicates that 6 c (-?Vhi) measured from the smoothed 
b distribution can be considered as a lower limit on the real 
b c (Nni), while b c (Nn 1 ) from the iterative power law fit can 
be considered as an upper limit on the real b c (Nn 1 ). 
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Table 4. The power law fit to the iVm-fr distributions 



<z> 






Sample A 








Sample B 






log(c , p ) 


(T T - 1) 


&c,i3.s (km s" 1 ) 10 


# a 


log(co.p) 


(It - 1) 


&o,i3.5 (km s" 1 ) 15 


1.61 


47 


-0.92 ±0.09 


0.17 ±0.01 


24.5 ± 2.4 


56 


-0.87 ±0.10 


0.17 ±0.01 


24.9 ± 1.1 


1.98 


103 


-0.49 ±0.11 


0.14 ±0.01 


21.4 ±0.5 


146 


-0.43 ± 0.12 


0.13 ±0.01 


21.4 ±0.6 


2.13 


139 


0.11 ±0.08 


0.09 ±0.01 


20.1 ±0.6 


181 


-0.70 ±0.11 


0.15 ±0.01 


20.6 ±0.5 


2.66 


140 


-0.55 ±0.10 


0.14 ±0.01 


19.6 ± 1.1 


204 


-0.13 ±0.10 


0.10 ±0.01 


18.7 ±0.5 


2.87 










223 


-0.89 ±0.08 


0.16 ±0.01 


20.4 ± 0.8 


3.75 


209 


0.51 ±0.06 


0.06 ±0.01 


18.7 ± 1.2 


271 


0.71 ±0.10 


0.04 ±0.01 


17.1 ±0.9 


2.1 c 


242 


-0.22 ± 0.08 


0.11 ±0.01 


19.6 ±0.6 


327 


-0.73 ±0.08 


0.15 ±0.01 


20.4 ± 0.5 


2.1 d 


156 


-0.64 ±0.13 


0.15 ±0.01 


21.2 ±0.7 


187 


-0.78 ±0.12 


0.16 ±0.01 


20.9 ±0.6 


3.75 c 


188 


0.09 ±0.10 


0.09 ±0.01 


18.9 ±1.0 


233 


-0.06 ±0.09 


0.10 ±0.01 


18.1 ±0.9 



The number of absorption lines used for the fit. 

The errors were estimated from the difference between the mean 6 c ,i3.5 and the minimum/maximum 6 c ,i3.5 from the 200 
bootstrap realizations for the given Nni~b pairs. The typical l a of t he & c ,i3.5 distribution from the 200 bootstrap realizations 
is ~ 0.3 km sec -1 , which is underestimated (cf. Schaye et al. 2000). 

Combined line lists from J2233-606 and HE2218-2817. The power law fit was carried out for AT H , = lO 12 5 " 14 ' 5 cm -2 . 
Combined line lists from J2233-606 and HE2218-2817. The power law fit was carried out for 7V Hl = lO 13 " 145 cm" 2 . 
For N Hl = 10 13-14 5 cm" 2 . 



Table 


5. The power 


law fit to the smoothed 


tions 








<Z> 


l°g(co, s ) 


(It - 1) 


fec,13.5 








(km s" 1 ) 


1.61 


-0.92 ±0.13 


0.16 ±0.04 


20.3 ± 1.1 


1.98 


-1.42 ±0.02 


0.20 ±0.01 


19.1 ± 1.0 


2.13 


-0.93 ±0.11 


0.16 ±0.03 


18.8 ± 1.1 


2.66 


-0.73 ±0.22 


0.14 ±0.06 


16.6 ± 1.2 


2.87 a 


-1.42 ±0.08 


0.20 ±0.02 


19.0 ± 1.1 


3.75 


0.16 ±0.07 


0.08 ±0.02 


16.4 ± 1.1 


2.1 b 


-0.71 ±0.14 


0.15 ±0.04 


19.2 ± 1.0 


2.1 c 


-0.50 ±0.34 


0.13 ±0.09 


19.7 ± 1.1 


3.75 d 


-0.30 ±0.13 


0.11 ±0.03 


15.8 ± 1.0 



For Sample B since Sample A cannot be defined. 



b At Nm = 10 1 ' 

2817. 
c At Nm = 10 13 " 
d At 7V HI = 10 13 ^ 



cm"' for J2233-606 and HE2218- 



cm" 2 for J2233-606 and HE2218-2817. 



As with the iterative power law fit, 6 c ,i3.5 measured 
from the smoothed b distribution increases continuously 
as z decreases, except at < z > = 2.87, where 6 c ,i3.5 is 
higher than at the adjacent redshifts (see Sect. 6.2 for 
further discussion). The slope (IV — 1) measured from 
the smoothed b distribution also does not show any well- 
defined trend with z. 

4.3.4. The b distributions 

Assuming that absorption lines arise from local optical 
depth (r) peaks and that In r is a Gaussian random vari- 



able, Hui & Rutledge ( 1999 ) derived a single-parameter b 
distribution: 



dn j db 



Bhr -jp- exp( 



b 4 



), 



(3) 



where n is the number of absorption lines, -Bhr is a con- 
stant and b a is a parameter determined by the average 
amplitude of the fluctuations and the effective smoothing 
scale. 

Fig. [n^ shows the observed b distributions at each z. 
The noticeable difference between Sample A (solid lines) 
and Sample B (dot-dashed lines) occurs at b < 20 km s" 1 
or b > 40 km s . These unphysical lines are usually in- 
troduced by VPFIT to fit noises so that the overall profile 
of Hi forest complexes could be improved. The dashed 
line represents the best-fitting Hui-Rutledge b distribu- 
tion, while the dotted line represents the b parameter for 
which the Hui-Rutledge b distribution function vanishes to 
10~ 4 , &hr, i-e. the truncated b value for the Hui-Rutledge 
b distribution function. The parameter &hr cannot be con- 
sidered equivalent to the cutoff b c since it is derived from 
the b distribution without assuming the b c dependence on 
TVhi. It is more sensitive to smaller b values in the b dis- 
tribution, which are in general coupled with lower iVni- 
Table || lists the relevant parameters describing the Hui- 
Rutledge b distribution for Sample A, such as the constant 
£?hr, b a , &hr and the median b values at different column 
density ranges. 

It is hard to specify subtle differences among the b 
distributions: while the modal b value and the 6hr value 
have a minimum at < z > — 3.75, they have a maximum 
at < z > — 2.87. The < z > = 2.87 forest also has the 
broadest b distribution. However, this large a(b) could be 
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Fig. 8. The logiVni-fc diagrams with the power law fits 
at z ~ 2.1 (only for J2233-606 and HE2217-2818; small 
filled circles) and at <z> = 3.75 (open squares). At both 
redshifts, the data points with errors less than 25% in both 
Nh i and b are displayed. The upper panel shows the fit to 
the column density range iVm = 10 12 - 5 - 14 - 5 rm -2 , while 
the lower panel is to the range Njii = 10 13 ~ 14 - 5 rar 2 . The 
vertical dotted lines in each panel represent the column 
density range over which the power law fit was carried 
out. The solid line shows the power law fit at z ~ 2.1, 
and the dashed line at z — 3.75. A deficit of lines with 
Nm < 10 13 cm~ 2 and b < 20 km s _1 observed at <z> = 
3.75 is due in part to the fact that all the fitted lines with 
./Vhi < fO 13 cm -2 and b < 20 km s _1 have errors larger 
than 25 % and in part to the severe line blending which 
limits the detection of weak lines. 
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Fig. 9. The log Nm-b diagrams for Sample A at each 
z (Sample B at < z > = 2.87) with the robust least- 
squares power law fit to the smoothed b distribution with 
the Gaussian smoothing constant 3 km s . Small open 
squares represent the data points for Sample A (Sample 
B at <z> = 2.87), while filled circles represent the cutoff 
b values estimated from the smoothed b distributions with 
~ 30 lines in each group (see the text for the details). The 
solid line represents the robust least-squares power law fit 
to filled circles. 

Table 6. The parameters of the b distributions 



<z> 


Bub 




6hr 


°med 


°med 




(km s- 1 ) 


(km s" 1 ) 


(km s- 1 ) 


(km s- 1 ) 


1.61 


7.37 


23.01 


12.61 


28.14 


34.56 


1.98 


7.98 


23.83 


13.04 


26.04 


29.10 


2.13 


7.46 


23.61 


12.95 


25.34 


29.57 


2.66 


6.82 


24.09 


13.25 


28.30 


30.10 


2.87 


7.09 


27.75 


15.30 


28.74 


34.05 


3.75 


6.72 


22.41 


12.31 


28.90 


30.70 


a For lines with iVm = 10 13 


x " 14 cm" 2 


from Sample 


A (Sam- 


pie 


B at <2> 


= 2.87). 








b For lines with AT Hl = 10 13 


8 " 16 cm" 2 


from Sample 


A (Sam- 


pie 


B at < z > 


= 2.87). 









in part due to a different fitting program and in part due 
to a lack of information on the errors. Other parameters, 
such as b(j , &hr, and cr(b), appear to be fairly constant 
with z. 

Fig. pd| shows the b distribution with z. This diagram 
does not assume a b c dependence on iVm, but is sensitive 
to a local 6 c (Ahi) variance. At z < 3.1, there is no clear 
indication of the behavior of the lower cutoff b values as 
a function of z. However, there is a clear indication of a 
trend with z of the lower cutoff b values over 3.5 < z < 3.9. 
In Fig. [ll], there are distinct regions at 1.8 < z < 2.4. 



The apparent cutoff values in&at2.2<z<2.4 and at 
1.8 < z < 1.9 are clearly higher than at 1.9 < z < 2.2. The 
2.2 < z < 2.4 region towards HE2217-2818 corresponds 
to a - 44 h,- 1 Mpc void (the region B in Fig. |l|), which 
might suggest enhanced ionization due to a nearby QSO 



or processes of galaxy formation (Theuns et al. 2000a ) 



4-4- The two-point velocity correlation function 

The Lya forest contains information on the large-scale 
matter distribution and the simplest way to study it is 
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Fig. 10. The 6 distribution of the Lya forest at each z. 
While solid lines are for Sample A, dot-dashed lines are for 
Sample B (no Sample A at < z > = 3.75). The dashed line 
and the dotted line represent the best-fitting Hui-Rutledge 
function and &hr, respectively. The 6 mo d value in the pan- 
els indicates the median b value at the corresponding z for 
Sample A. The number in parentheses indicates the lcr(6) 
value from the Gaussian distribution. 



to compute the two-point velocity correlation function, 
£(Au). The correlation function compares the observed 
number of pairs (iV b s ) with the expected number of 
pairs (A exp ) from a random distribution in a given ve- 
locity bin (Av): £(A«) = N ohs (Av)/N cxp (Av) ~ 1, where 
Av — c (z 2 — zi)/[l + (z2 + zi)/2], z\ and z 2 are redsh ifts of 
two lines and c is the speed of light (Cristiani et al. 1995; 
Cristiani et al. |997} Kim et al. |l997j ). 

Studies of the correlation function of the Lya forest 
have generally led to conflicting results even at simil ar z 
Some studies find a lack of clus terin g (Sargent et al 
at 1.7 < z < 3.3; Rauch et al 
al. 1994 at z ~ 4), while othe rs hnd clustering at scales 
Av < 350 km s" 1 (Webb |l987| at 1.9 < z < 2.8; Hu et al. 
19951 at z ~ 2.8; Kulkarni et al. Jl996| at z ~ 1.9; Lu et al. 
199q at z ~ 3.7; Cristiani et al. |1997| at z ~ 3.3). 
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Fig. 11. The b distribution as a function of z for Sample 
A (Sample B at <z> — 3.75). The horizontal dashed line 
indicates 6 = 20 k m s -1 which is a Ah r independent b c at 
z ~ 2.8 (Hu et al. 1995 ). Circles, crosses, diamonds, stars, 



triangles and squares are from HE0515-4414, J2233-606, 
HE2217-2818, HS1946+7658, Q0302-003 and Q0000-263, 
respectively. There is an indication of increasing b c with 
decreasing z at z ~ 3.7. At z < 3.1, b c is not clearly 
defined. 



(HE0515-4414, J2233-606, and HE2217-2818), 2.5 < z < 
3.1 (HS1946+7658 and Q0302-003) and 3.5 < z < 4.0 
(Q0000-263). 

In our approach A^xp was estimated averaging 1 000 
numerical simulations of the observed number of lines, 
trying to account for relevant cosmological and observa- 
tional effects. In particular a set of lines was randomly 
generated in the same redshift interval as the data ac- 
cording to the cosmological distribution oc (1 + z) 7 , with 
7 = 2.4 (see Sect. 4.2). The results are not sensitive to 
the value of 7 adopted and even a flat distribution (i.e. 
7 = 0) gives values of £ that differ typically by less than 
0.02. Line blanketing of weak lines due to strong complexes 
was also accounted for. Lines with too small velocity split- 
tings, compared with the finite resolution or the intrinsic 
blending due to the typical line widths-the so-called "line- 



blanketing" effect (Giallongo et al. 1996), were excluded 
in the estimates of N exp . 

Clustering is clearly detected at low redshift: at 1.5 < 
z < 2.4 in the 100 km s _1 bin, we measure £ = 0.4 ± 



1992 at z - 3; Williger et 0.1 for lines with iV H i k, 



of increasing amplitude with increasing column density: 
in the same redshift range £ = 0.35 ± 0.08 for lines with 
iVm k 10 12 5 cm -2 . The trend is not significant but agrees 
with the behavior observed at higher redshifts (Cristiani 



Fig. p. 2| shows the velocity correlation strength at Av < 
4 000 km s^ 1 . To obtain sufficient statistics, the analy- 
sis was carried out in three redshift bins: 1.5 < z < 2.4 



et al. 1997; Kim et al. 1997). Unfortunately the number 



of lines observed in the interval 1.5 < z < 2.4 does not 
allow us to extend the analysis to higher column densities, 



Tae-Sun Kim, Stefano Cristiani & Sandro D'Odorico: The Lya forest at 1.5 < z < 4' 



15 



0.4 
0.2 


-0.2 

0.4 
0.2 

> 

o 

-0.2 

0.4 
0.2 


-0.2 



— , I | I I I i | I I I I | ,— 

log(N HI ) > 12.7 1.5 < z < 2.4 



. i j. 



i _£i 



H 1 1 h 



H 1 h 



H 1 1 h 



H 1 h 



2.5 < z < 3.1 



XL X*X~_ _ ,i A _ V *LTT» _ _A ^ _ _ _ /TVT IK _ _ 

:x_ LLrxxxit!!: ix zx_x_ r_~_z.~ ~ ~ V 



H 1 h 



3.5 < z < 4.0 



Y • — — - 



HE0515-4414 



1000 2000 3000 

Av (km s _1 ) 

Fig. 12. Evolution of the two-point correlation function 
with redshift for Lya lines with column densities above 
TVhi = lO 12 7 cm~ 2 . Short-dashed and long-dashed lines 
represent the la and 2a Poisson errors. 



although groups of strong lines are occasionally evident 
(e.g. the range 3230-3270 A in HE0515-4414). 

The amplitude of the correlation at 100 km s _1 de- 
creases significantly with increasing redshift from 0.4 ±0.1 
at 1.5 < z < 2.4, to 0.14 ± 0.06 at 2.5 < z < 3.1 and 
0.09 ±0.07 at 3.5 < z < 4.0. 



4-5. Voids 



Voids along the three low-redshift lines of sight were 
searched for. For comparison with previous results (Car- 



swell & Rees |1987[ Crotts |1987| ; Ostriker ct al. |198S| ), we 
identify a void as a region without any absorption stronger 
than AThi ~ 10 13 5 cm~ 2 over a comoving size of at least 
30 h^ 1 Mpc (assuming q — 0.1). 

Figs. [l3|-|l4] show the voids detected in the spectrum of 
HE0515-4414 and HE2217-2818, respectively. No signifi- 
cant void was found in the spectrum of J2233-606. The 
wavelength range used for searching for voids has been 
selected to be redward of the QSO's Ly/3 emission line 
and 3 000 km s _1 blueward of the QSO's Lya emission to 
avoid the proximity effect. The wavelength range searched 




3100 3150 3200 3250 
Vacuum Wavelength (A) 

Fig. 13. The spectrum of HE0515-4414 with the void at 
z = 1.570. See the text for the details. 

for voids is larger than that used to study the Lya forest 
in other sections. 

Table lists the dimensions of the voids, as well as the 
probability of finding a void larger than their comoving 
size. The probability was calculated assuming a Poisson 
distribution of the local forest. In this case, the probability 
of finding a void larger than a given size x gap is P> (x gap ) — 
1 — (1 — exp _:Csap )™, where x gap is the line interval in the 
unit of the local mean line interval and n is the number of 
lines with AT Hl > 10 13 5 cm" 2 (Ostriker et al. |198SQ . The 
joint probability of finding two voids with a size larger 
than 40/iX 1 Mpc at z ~ 2, as observed in the spectrum 
of HE2217-2818, is of the order of 2 x 10~ 4 . The results 
correspond very well to the probability estimates derived 
from the simulations described above. 

There are different ways to produce a void in the for- 
est: a large fluctuation in the gas density of absorbers, en- 
hanced UV ionizing radiation from nearby faint QSOs or 
star- forming galaxies, feedback processes (including shock 
heati ng) from gal axy f ormation (Dobr zycki fc Bechtold 
199lfc Heap et al. |2000| ; Theuns et al. |2000a| ). We recall 



here that the void B in the spectrum of HE2217-2818 cor- 
responds to a region of above-average Doppler parameter 
(Sect. 4.3.3). It will be interesting to carry out deep imag- 
ing around HE22 17-28 18 to identify QSO candidates and 
investigate whether a local ionizing source is responsible 
for the ~ 50 h~ x Mpc voids. 

5. The flux statistics on the Lya forest 

The traditional Voigt profile fitting analysis is limited by 
two major drawbacks. First, there is no unique solution. 
Although the \ 2 minimization is normally applied to the 
fit, different fitting programs produce slightly different re- 
sults. Even using the same program, different x 2 thresh- 
olds lead to different numbers of lines when line blending 
is severe. As the resolution and S/N increase, many for- 
est lines show various degrees of departure from the Voigt 
profile. This departure can be fitted by adding physically 
improbable narrow components in high S/N data to im- 
prove the overall fit, while the same profile can usually be 
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Table 7. Voids at z ~ 2 



QSO 




Region 


Wavelength 


<z> 


Az 


Comoving size a 


Kgap 


P>(a;ga P ) 








(A) 






(ft" 1 Mpc) 






HE0515- 


-4414 


A 


3088-3161 


1.570 


0.060 


61.1 


5.7 


0.045 


HE2217- 


-2818 


A 


3504-3579 


1.913 


0.062 


54.3 


8.4 


0.012 


HE2217- 


-2818 


B 


3878-3946 


2.218 


0.056 


43.5 


8.0 


0.018 



For q = 0.1. 



03 




3500 3600 3700 3800 3900 4000 4100 
Vacuum Wavelength (A) 

Fig. 14. The spectrum of HE2217-2818 with the void re- 
gions. The voids are indicated as A at z = 1.913 and B at 
z = 2.218. See the text for the details. 



fitted by one broad component in low S/N data (Schaye 
et al. |1999| ; Theuns & Zaroubi |2000[ ). Thus, it becomes 
harder to compare fitted line parameters from different 
observations and different fitting programs, as the resolu- 
tion and S/N increase. 

Second, there is no guarantee that absorption lines 
have a Voigt profile. Although searches for characteris- 
tic non-Voigt signatures have been inconclusive (Rauch 
1996); Outram et al. 1999b| ), numerical simulations pre- 



dict non-Voigt profiles of absorption lines (Outram et al. 
200C). In numerical simulations, lines are broadened by a 
combination of gas temperature, peculiar velocity, Hub- 
ble expansion, and Jeans smoothing (Bryan et al 



1999 



Theuns & Zaroubi 2000| ). In conclusion, the profile fitting 
method should be regarded as a powerful parameteriza- 
tion of the spectrum, not a representation of the physics 



governing forest clouds (Rauch et al. 1992; Hu et al. 1995 
Kirkman & Tytler |1997|; Bryan et al. |g9§ Schaye et al. 



1999; Theuns & Zaroubi 2000) 



In order to avoid the non-uniqueness of profile fitting 
analyses and to allow more straightforward comparisons 
with theoretical predictions, a direct use of observed spec- 
tra of the Lya forest has been explored. The most straight- 
forward way to characterize the observed spectra is to use 
the N-point functions of the transmitted flux, F, or the ob- 



et al 
et al 



1997 



2000 



serve d opt ical depth, r (Mira lda-Escude e t al. 1997 ; Rauch 
Zhang et al. |l998|; B ryan et al. |1999[ Machacek 
Theuns et al. |2000bp . In the following, we con- 



sider the one-point function and the two-point function of 
the transmitted flux as well as other statistical measures 
such as the line count and the optical depth correlation 
function. 

For comparison with the forest at z > 2.4, we gener- 
ated artificial spectra from the published line lists. These 
lists are the same as those used in the Voigt profile fit- 
ting analysis in Sect. 4 (see Table |J). We added Gaussian 
noise to the artificial spectra, according to the quoted S/N. 
However, as Theuns et al. ( [2000b ) stated, Gaussian noise 
independent of flux does not represent the real, observed 
S/N. The difference becomes more evident at F ~ (satu- 
rated regions), where Gaussian noise produces larger fluc- 
tuations than what is observed. Unfortunately, without 
details on the number of individual spectra for a given 
wavelength range, data reduction and the normalization 
of the spectra, the noise at the bottom of saturated lines 
cannot be simulated correctly just from the published S/N 
and the given CCD readout noise. 

We tested our approach with two extreme cases: ar- 
tificial spectra without noise and artificial spectra with 
Gaussian noise. Except for F~0at<z> = 3.75, the 
results from these two cases do not differ significantly, in 
particular, at 0.1 < F < 0.8 where most results are con- 
sidered. We use the results from the spectra generated 
with Gaussian noise in this study. However, we mention 
any difference between these two extreme cases when it 
becomes noticeable. Also note that most simulations do 
not include noise in their analysis and, in a sense noise- 
free spectra would be more appropriate to be compared 
with the results from simulations. Due to the imperfect 
simulation of the noise in the artificially generated spec- 
tra, the results drawn from the flux statistics should be 
taken qualitatively at F ~ and F ~ 1. 

In order to be consistent with our profile fitting analy- 
sis, we excluded high-column-density systems with Vjji > 
10 16 cm -2 from the spectra. Note that including these 
high-column-density systems does not change the general 
qualitative conclusions from the flux statistics. It only 
changes the quantitative results considerably at z ~ 3.7 
when F becomes close to due to the higher number of 
pixels with F ~ at z ~ 3.7. However, for the estimation 
of the mean Hi opacity, we used the whole regions in- 
cluding high-column-density absorption systems in order 
to compare with the previous results from the literature 
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(Press et al. 1993; Rauch et al. 1997). We also remind the 
reader to be particularly cautious in using the line list of 
Q0302-003 at z ~ 2.9, which is the only one not generated 
by VPFIT. As for the Voigt fitting analysis, a systematic 
difference in iVni and b could change the results in this 
Sect, at z ~ 2.9. 



5.1. The mean Hi opacity 

The Hi opacity, rm(A), is defined as Tm(A) = 
— ln(F\/F c ), where A is the observed wavelength, Fa is 
the observed flux at A, and F c is the continuum flux at A. 
Since the opacity scales logarithmically, the mean opac- 
ity cannot be measured accurately when F\ « 0. The 
effective opacity, r e g, is typically used in place of thi: 
exp _Tcff = <exp _r >, where < > indicates the mean value 
averaged over wavelength. Note that the mean Hi opacity 
and the effective Hi opacity are different quantities. How- 
ever, we refer to the estimated r e gf values as the "mean 
H I opacity" in this study. 

Fig. [l5| shows our fm measures (in fact, the effec- 
tive optical depth T e ff), together with other opacity mea- 
sures compiled from the literature. Filled circles are the 
thi measures from the UVES data including high-column- 
density absorption systems. Note that all the other opac- 
ities from the literature (except the Press et al. ( 1993 ) 
thi measurements) were measured including high-column- 
density regions in the spectra except damped Lya sys- 
tems. The Press et al. measure (dotted line) was derived 
at low resolution, including damped Lya systems. Our 
experiments with the UVES data show that there is no 
noticeable difference between the thi measured from the 
spectra generated from the line lists and the ?m from 
the observed spectra. This indicates that very weak lines 
do not significantly contribute to rm- Therefore, the esti- 
mated th i from the spectra generated from the published 
line lists using high resolution, high S/N data can be con- 
sidered to be reliable. In Table g, we list the estimated thi 
values when these values are not given in numeric form in 
the relevant references of Fig. U3. 



The opacity measures are dependent on the continuum 
fitting. For low-resolution data, the continuum is usually 
extrapolated from longward of the QSO's Lya emission 
line, typically resulting in an overestimation of the con- 
tinuum, i.e. an overestimated fm. On the other hand, 
the local continuum fitting generally adopted for high- 
resolution data may result in an underestimation of the 
continuum, i.e. an underestimated thi- Therefore, it is not 
surprising that the thi (dotted line) from low-resolution 
data by Press et al. (1993) is higher than any other mea- 
surements. Small triangles from Zuo & Lu (1993) were 



estimated from the published line lists using intermediate 
resolution data, which usually do not include low column 
density lines. Thus, the Zuo & Lu estimates are about a 
factor of 2 lower than the PRS formula at z ~ 3. Other 




Fig. 15. The Hi opacity as a function of (1 + z). Filled cir- 
cles represent the mean H I opacity from the UVES data. 
Sma ll tria ngles and diamonds represent thi from Zuo & 
Lu ( |1993| ; adding up the equivalent width s from the pub- 
lished line lists) and Rauch et al. ( 1997 ; we used their 
t|™' uncorrected for continuum fitting uncertainty since 
all other observations were not corrected for continuum 
fitting uncertainties), respectively. Open circles, the star, 
the large square, the small diamond and the small square 
are estimated from the spectra generate d art ificially us- 
ing th e line list provid ed by Hu et al. ( 1995 ), Kulkarni 
et al. (|1996D , Lu et al. (|1996|) , Kirkman & Tytler ([L997D , 
and Outram et al. ( 1999a ), res pectiv ely. Large triangles 
were read from Schaye et al. ( |2000 ). All the x-axis er- 
ror bars represent the z range used for the ?m estimates, 
while the y-axis error bars were estimated from simply 
changing the adopted continuum by ±5%, except for large 
triangles whose y-axis errors were read from Schaye et 
al. ( |2000D . The cross at z = 4.7 (Songaila et al. 
and the thick diamond at z — 5.34 (Spinrad et al 



1999) 



1998) 



are from the Keckll/LRIS data. The dotted line repre- 
sents the commonly used formula by Press et al. (1993): 
thi(z) = 0.0037(1 + z) 3 - 46 , given by the conventional 
power law fit. The solid line represents the power law fit 
to the UVES and the HIRES data. 

observations fall inbetween the PRS formula and the Zuo 
& Lu values. 

There is a scatter in thi at similar z, even though 
opacities are estimated from data with the same instru- 
ment configuration and reduction, such as two filled circles 
at z ~ 2.1 from UVES and three open circles at z ~ 2.8 
from HIRES. This scatter could result from measurement 
errors, such as the continuum fitting, and from a cosmic 
spatial variance. The continuum fitting becomes unreliable 
at z > 3 due to severe line blending. On the other hand, 
the cosmic variance is present at all z. For example, the 
difference between two UVES measurements at z ~ 2.1 
(filled circles) is due to the two high-column-density sys- 
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Table 8. The mean Hi opacity 



QSO 



<z> 



Tin 



Ref. 



0.086% J49 51 
0.064°_g 4 9 51 

0.1128:8?? 

0.161°_-8 4 9 51 
0.142°_'8 4 o 9 5i 

n i qi 0.049 
U.ldl_o.051 

0.234°_-8 4 o 9 5i 
0.298°_-8 4 o 9 5i 
0.275°_-8 4 o 9 5i 
0.289°_-8 4 o 9 5i 

1.049 
-0.051 
;0.04 
-0.04 

0.83°_8 6 08 



HE0515-4414 

Q1331+170 

Q 1100-264 

J2233-606 

J2233-606 

HE2217-2818 

HS1946+7658 

Q0636+680 

Q0302-003 

Q0014+813 

Q0000-263 

Q2237-061 

Q2237-061 



1.61 
1.85 
1.96 
1.98 
2.04 
2.13 
2.66 
2.80 
2.87 
2.97 
3.75 
3.84 
4.31 



1.54-1.68 

1.68- 2.01 
1.85-2.09 
1.80-2.17 
1.92-2.17 
1.92-2.33 
2.50-2.81 
2.58-3.02 
2.63-3.11 
2.74-3.20 
3.48-4.02 

3.69- 4.02 
4.15-4.43 



0.75° 



This study 

Kulkarni et al. (|l996|) 
Schaye et al. ( EOOOD 
This stud y 
Outram et al. ( 1999ej ) 
This stu dy 
Kirkman & Ty tler (|l997| ), HIRES 
Hu et al. (11995), HIRES 
Hu et al. (1995), HIRES 
Hu et al. (1995), HIRES 
Lu et al. p996|) HIRES 
Schaye et al. (pOOOl), HIRES 
Schaye et al. (2000), HIRES 



terns with iVm > 10 17 cm 2 towards J2233-606. When 5.2. The one-point function of the flux 



these high-column-density systems are excluded, the opac- 
ity at < z > = 1.98 (J2233-606) becomes 0.123, which is 
similar to the one at < z > = 2.13 (HE2217-2818). Simi- 
larly, the higher opacity at z ~ 2.80 towards Q0636+680 
is due to several high-column-density clouds on this line 
of sight, while two other lines of sight at z ~ 2.8 are rel- 
atively devoid of high-column-density clouds. Despite the 
same HIRES configurations, the thi values from Hu et al. 
Q1995Q and from Kirkman & Tytler ( |1997| ) are a factor of 
1.2 lower than the Rauch et al. values. Since the line lists 
for the Rauch et al. QSO sample are not published, we 
cannot test whether the presence of high-column-density 
systems in their sample causes the higher rm- 

The solid line represents the least-squares fit to the 
UVES and the HIRES data: r Hl (z) = 0.0030 ±0.0008 (1 + 
z) 3 - 43±017 . The new UVES data at 1.5 < z < 2.4 suggest 
that thi can be well approximated by a single power law 
at 1.5 < z < 4. Two thi measures at z > 4.5 (cross 
and thick diamond) from the Keckll/LRIS data suggest 
that thi might be significantly higher than extrapolated 
from z < 4. However, these values were derived from low- 
resolution data, which usually overestimate ?m. In fact, 
they correspond better to the Press et al. formula, which 
was also derived from low-resolution data. Without more 
high-resolution data at higher z, it is premature to con- 
clude that the fm evolution at z > 4 departs significantly 



from a single power law (cf . Schaye et al. 200C ) . 

In the standard Friedmann universe with the cos- 
mological constant A = 0, if the baryon overdensity 5 
is 6 < 10, the mean Hi opacity can be expressed by 



thi oc (1 + zy- 6 /Tm(z) (Machaceket al. 2000). Although 



there are a few high-column-density systems in Fig. |15|, 
our result on r H i, t H i oc (1 + z ) 3 - 43 ± 017 a t 1.5 < z < 4, 
is in good agreement with the predicted power law index. 
This suggests that Thi(^) does not strongly evolve over 
this z range. 



The one-point function of the flux (or the probability den- 
sity distribution function of the transmitted flux), P(F), 
is simply the number of pixels which have a flux between 
F and F + dF for a given flux F over the entire number of 
pixels per dF. In other words, it is the probability density 
to find a pixel at a given F (Miralda-Escude et al. [1997 



Rauch et al. |1997|; Br yan et al. [1999|: Machacek et al. |2000 
Theuns et al. |2000b[ ). 

Fig. |l| shows P(F) as a function of F. The one-point 
function of the flux at F < and F > 1 from observations 
is due to observational and continuum fitting uncertain- 
ties. The non-smooth P(F) at < z > = 1.61 is due to 
the small number of pixels used to calculate P(F). The 
wider P(F) profiles at F ~ and at z > 2.4 compared 
with at z < 2.4 are due to the characteristics of Gaussian 



noise in the spectra generated (Theuns et al. 2000b ). For 
the spectra generated without noise, the P(F) profiles at 
z > 2.4 are narrower with higher amplitude at F ~ 0, but 
do not differ significantly from the P(F) profiles in Fig. |l^ 
at 0.2 < F < 0.8. 

The flattening towards F^lat<z> = 3.75 is due 
to the smaller number of F ~ 1 pixels from the increasing 
forest number density at higher z. The one-point functions 
of the flux at z ~ 2.1 (for J2233-606 and HE2217-2818) 
are a factor of 1.3 hig her th an the sCDM model simulated 
by Machacek et al. fl2000|) at 0.2 < F < 0.6 (this lower 
simulated P(F) is in agreement with their lower mean Hi 
opacity at z ~ 2). 

After smoothing P(F) at < z > = 1.61 over a 

dF = 0.14 bin, P(F) at F = 0.2 becomes P(F) F=0 .2 oc 
(1 + z )3.86±o.54^ Thc flux F = 2 corresponds t0 A r Hi _ 

3.2 x 10 13 cm~ 2 (in this study, we assume b = 30 km 
s- 1 ). At F — 0.6 (JV H , = 1.0 x 10 13 cm" 2 ), P(F) F=0 . 6 oc 
(1 + z) 3 " 26±0 ' 21 . In short, the probability density of finding 
strong absorption lines (smaller F) shows a steeper slope 
than that of finding weak absorption lines (larger F) . This 
is in good agreement with the result from the Voigt profile 
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Fig. 16. The one-point function of the flux, P(F), as 
a function of F at different z. Strong absorption lines 
(smaller F) disappear more rapidly than weak absorption 
lines (larger F) as z decreases. 

fitting analysis: the higher column density forest disap- 
pears more rapidly than the lower column density forest 
as z decreases. 



5. 3. The two-point function of the flux 

The two-point function of the flux, P{F\, F2, Av), is the 
probability of two pixels with the Av velocity separation 
having normalized fluxes F\ and F2 . It is usually expressed 
as 

AF(Aw,<5Fi) 

= I ! {F 1 -F 2 )P(F 1 ,F 2 ,Av)dF 2 

JSFx YJ-00 



dF 1 /8F l , (4) 



where AF(Av,SFi) is the mean flux difference between 
two pixels with F\ and F 2 , which are separated by Av. 
At Av > 200 km s _1 , ~AF(Av,6F 1 ) is just the dif- 
ference between the mean transmitted flux and F±. At 
30 < Av < 100 km s" 1 , the shape of AF(Au, 6 Fx) con- 
tains information on the profile shape of absorption lines. 
However, great care has to be exercised in the interpreta- 
tion since AF(Av,SF±) is an average quantity (Miralda- 
Escu de et al. |1997[ ; Bryan et al. [199$ Machacek et al. 
2000| ; Thcuns et al. |2000b[ ). 

Fig. shows AF(Av, &F\) as a function of Av at dif- 
ferent z. The left-hand panel shows the results for the 
observed spectra at z < 2.4 and the spectra generated 
with noise at z > 2.4, while the right-hand panel shows 
the results for the spectra generated from the line lists 
for all z without noise. There is no noticeable difference 
between the two panels. To b e com parable with the sim- 
ulations by Machacek et al. ( 2000| ), the flux range was 
chosen to be —0.1 < F < 0.1, which corresponds to 




10 100 
Ai> (km s~ 1 ) 



10 100 
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Fig. 17. The mean flux difference, AF(Av,SF 1 ), as a 
function of Av at different z. Since AF(Av, 6F\) is an av- 
eraged quantity, it is not obvious to interpret AF(Av, SFi) 
to be a strong indicator of the temperature of the IGM. 



JVhi > 10 13 ' 77 cm -2 . Over this column density range, ab- 
sorption lines are in general saturated and belong to a H I 
complex, rather than being isolated. 

The overall z-dependence on AF(Av, SFi) is not clear. 
At Av < 30 km s _1 , AF(Av,8Fi) shows almost identi- 
cal profiles at all z except at < z > — 1.61 and 2.87. At 
< z > = 1.61 and 2.87, A~F(Av,6Fi) is wider than at 
any other z. At Av > 30 km s _1 , there is a tendency 
for the width of AF(Av,SFi) to become narrower as z 



decreases. Following Machacek et al. ( 2000 ), we measured 
Au(0.3), the width of AF(Av, SFi) at which 7AT(Aw, SFi) 
becomes 0.3. At < z > = 1.61, 1.98, 2.13, 2.66, 2.87, and 
3.75, A«(0.3) is 51.06 (51.61), 44.25 (43.15), 41.62 (44.27), 
47.05 (48.54), 49.99 (50.79) and 50.95 (52.11) km s" 1 (the 
number in parentheses is for the spectrum generated with- 
out noise), respectively. The observed Av(0.3) values are 
larger than the predicted Au(0.3) values, ~ 35 km s -1 . 



from the models considered by Machacek et al. (200C) at 
z~3. 

A weak trend of decreasing Au(0.3) with decreasing 
z might suggest that the gas temperature decreases as 
z decreases. In fact, Theuns et al. (2000b) note that at 
a given z a simulation with a hotter gas temperature 
shows a wider AF(Av, 6F±) profile than a simulation with 
a lower gas temperature. However, numerical simulations 
also show that higher b values at higher z may be a result 
of other physical processes. Variations of AF(Av, SFi) as 
a function of z can be a result of increasing line blending 
at higher z as well as a change in the gas temperature. 

5.4. Line counts of the Lya forest 

The line count at a given flux F t is defined a s the num- 
ber of regions below F t (Miralda-Escude et al. 1996 ). The 
number defined this way is more straightforward to deter- 
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Fig. 18. The line counts of the Lya forest. The upper 
panel shows the normalized flux as a function of the filling 
factor. The lower panel represents the line counts from the 
observed spectra and the spectra generated with noise, 
which were smoothed with a 20 km s^ 1 box-car function 
in order to decrease noise. 



mine than the conventional line counting from the profile 
fitting. 

The upper panel of Fig. |l8| shows the normalized flux 
of the observed spectra and the spectra generated with 
noise as a function of the filling factor. The filling fac- 
tor is the fraction of the spectrum occupied by the pixels 
whose normalized flux is smaller than a given F t . Since this 
method is sensitive to noise, the spectra were smoothed 
with a 20 km s _1 box-car function. Except at < z > = 3.75 
(the broad feature at the filling factor close to is due to 
the characteristics of Gaussian noise, i.e. larger root-mean- 
square fluctuations than the real, observed fluctuations at 
F ~ 0), the line counts are similar when the filling factor 
is 0.07-0.3. This range of the filling factor corresponds to 

0. 2 < F < 0.9 at z ~ 1.8 and to 0.05 < F < 0.8 at z ~ 2.7, 

1. e. pixels with a wide range of fluxes are considered. How- 
ever, the same filling factor range only probes < F < 0.3 
at < z > = 3.75, where only strong lines are counted. This 
results in the different line count at < z > — 3.75. When 
the filling factor is greater than 0.3, the line counts devi- 
ate from each other. This regime corresponds to F ~ 1, 
where noise distorts the true line counts. 

Fig. [l9] shows the line counts as a function of the fill- 
ing factor again. In this diagram, the line counts were 
calculated using the artificial spectra generated from the 
fitted line lists for each QSO without adding noise, i.e. 



they have an infinite S/N. Note that this process does 
not include weak lines, usually not present in the fitted 
line lists. Unlike the lower panel of Fig. [l8|, the curves de- 
scribing the line counts as a function of the filling factor 
show a similar progression as a function of z when a fill- 
ing factor is smaller than ~ 0.3, except for <z> — 2.87. 
The < z > — 2.87 forest also shows a slightly different 
behavior in the filling factor-flux diagram. The flux at a 
given filling factor increases continuously as z decreases 
except that the flux corresponding to a given filling fac- 
tor is larger at < z > — 2.87 than at < z > — 2.66 for a 
filling factor larger than 0.5. This might indicate the real 
cosmic variance in the structure of the Lya forest along 
the line of sight towards Q0302-003 (one known void to- 
wards Q0302-003 at z ~ 3.17 is not included in this study. 
Also note that the < z > = 2.13 forest includes one void 
region). However, a similar work done by Kim ( 1999] ) did 
not show any difference in the line counts as a function 
of z at 2.1 < z < 4 when the real observed spectra in- 
cluding Q0302-003 were used. It is highly unlikely that 
the QSO sample in Kim's work shows the same amount of 
systematic differences from the QSOs used in this study 
for all z. Therefore, we cast doubts on the cosmic vari- 
ance as a probable reason for the different line counts at 
< z > = 2.87. Rather, the line parameters of Q0302-003 
have not been obtained with the program VPFIT and this 
suggests that a different behavior of the Lya forest at 
z ~ 2.9 from the rest of the forest at different z, such as 
a higher 6 c ,i3.5 than at adjacent z, should be taken with 
caution. 

Since the filling factor is determined mainly by the 
Hubble expansion, the negligible ^-dependence of the line 
counts suggests that the evolution of the forest at 1.5 < 
z < 4 is drive n mai nly by the Hubble expansion (Miralda- 
Escude et al. 



1996) 



5. 5. The step optical depth correlation function 



Miralda-Escude et al. (1996) first introduced a correlation 
function using a pixel-by-pixel transmitted flux, which is 
more straightforward than the two-point velocity corre- 
lation function. Cen et al. (1998) developed this concept 
further. We analyzed the clustering properties of the Lya 



forest, following Cen et al.'s methods (1998). Among their 
newly defined correlation functions, we only consider the 
step optical depth correlation. In general, the trends we 
found from the step optical depth correlation function 
hold for the other correlation functions. 

The step optical depth correlation function £ T:S is de- 
fined as 



£t, s (Ai;) 



< r s (v + Av)t s (v) > 



L 



where the step optical depth is defined as 



t s (Av) = 



1, if Tobs (v)> 

0, if Tobs (v) < 



(5) 



(6) 
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Fig. 19. The line counts of the Lya forest. The upper 
panel shows the normalized flux as a function of the fill- 
ing factor. The lower panel represents the line counts of 
the spectra generated artificially from the line lists. The 
abrupt drops in the line counts at the larger filling factor 
are due to the fact that the large number of the F = 1 
pixels do not contribute to the line counts as the filling 
factor approaches 1. 



Fig. go] shows the step optical depth correlation func- 
tions £ T . s (Au) with r m ; n = 2 as a function of Aw. The 



2 corresponds to Am 



10 



13.6 



cm . Although 



it is difficult to detect any well-defined trend of the cor- 
relation strength at Av < 100 km s _1 in the two-point 
correlation functions, the step optical depth correlation 
functions show a strong clustering at v < 100 km s" 1 . 
Instead of the individual pixels, the two-point correlation 
function uses the fitted line lists from the Voigt profile 
fitting. The profile fitting usually does not deblend satu- 
rated lines (higher r and higher Am) and is less sensitive 
to fitting lines at small velocity separations. 

The step optical depth correlation strength increases 
as z decreases, except at < z > = 2.13 at Av < 50 km 
s -1 . The stronger correlation strength of the < z > = 2.13 
forest than the < z > = 1.98 forest at Aw < 30 km s _1 is 
in part caused by the 5 strong lines clustered at 3960-4005 
A. Also the < z > = 2.13 forest contains a void of 44 /i -1 
Mpc. 

The pixel-by-pixel correlation functions are sensitive to 
the S/N and resolution of the data. However, the general 
trend of increasing correlation strength with decreasing z 
holds at the different S/N and resolution from the exper- 
iments of degraded UVES spectra, although the degree of 
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Fig. 20. The step optical depth correlation functions. 
They are sensitive to the profile shapes of the forest at 
the smaller velocity separations for a given T m - ln . 

the correlation strengths gets weaker as the S/N and res- 
olution decrease. With higher T m j n values, the correlation 
strengths get stronger, but keep the same z-dependent 
trend. These results are in good agreement with previ- 
ous findings from the two-point velocity correlation func- 
tion: stronger forest lines are more strongly correlated and 
the correlation strengths increase with decreasing z over 
a given column density range (Cristiani et al. 1997; Kim 
et al. |1997|). 



6. Discussion 

6.1. The line number density per unit redshift of the Lya 
forest 

For a given column density range, combined HST and 
ground-based observations have provided evidence for a 
change in the line number density evolution at z ~ 1.7: a 
rapid evolution at z > 1.7 and a slow evolution at z < 1.7 



(Impey et al. 1996; Riediger et al. 1998. Weymann et al. 
|1998 ; Dave et al. 1999). These observations have led to a 
speculation of two distinct populations in the Lya forest: 
a rapidly evolving population which dominates at higher 
z and a slowly evolving population which dominates at 
lower z. 



Our results at Am 



10 



13.64-16 



cm 2 suggest that the 



transition from the stronger evolution to the weaker evo- 
lution in dn/dz occurs at z ~ 1.2 (UM18 fits in the pic- 
ture and is not an outlier), rather than at z ~ 1.7 as 
suggested by previous observations and numerical simula- 
tions. To be conservative, our results show that dn/dz at 
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iVm = 10 1364_16 cm~ 2 continues to decrease at a similar 
rate from z = 4 — > 1.5, with a suggestion of slowing down 
in the evolution towards lower z. 

The physics of the Lya forest at z > 2 is determined 
mainly by the Hubble expansion and the ionizing back- 
ground, J v (or the Hi photoionization rate, Thi)- If the 
forest is "fixed" in comoving coordinates for q n = 0.5 and 
A = 0, the observed number density of the Lya forest 
is proportional to [(1 + z) 5 T^] (z)]' 3 " 1 (1 + z)~ 3 / 2 , where 
[3 is from /(JVhi) oc iV«f . This implies that dn/dz for a 
given column density threshold decreases as z decreases 
(Miralda-Escude et al. |l996[ Dave et al. |l999| ). If we as- 
sume (3 ~ 1.5, dn/dz oc (1 + z)T^° 1 5 (z). For a constant 
Thi, this is much lower than the observed index of dn/dz, 
2.19±0.27, suggesting structure evolution and/or Lhi evo- 
lution in dn/dz. 

Recent numerical simulations suggest that a decrease 
of Thi at z < 2 plays a more important role to change 
the slope in dn/dz at z ~ 1.7 due to the decr easing QSO 



luminosity function at z < 2 (Dave et al. 1999; Riediger et 
al. |1998| ; Theuns et al. |199S| ; Zhang et al. |1998[ ). The dis- 
crepancy between our observations and simulations could 
be due to limited box sizes at z < 2 in most simulations 
(losing large-scale power), to numerical resolutions (un- 
derestimating thi or the number of lines at lower z) or 
to incorrect Thi- If we take the results from most simu- 
lations that Thi is the main drive of the slope change in 
dn/dz, this discrepancy could simply indicate that Tj^ in 
most simulations, i.e. the QSO-dominated Haardt-Madau 
(Haardt & Madau 1996| ), is underestimated at z < 2 
and that T^] at z < 2 changes more slowly than a QSO- 
dominated Tjjj, i.e. there is a non-negligible contribution 
from galaxies at z < 2. 



6.2. The z-evolution of b c (Nn 1 ) 

Assuming a truncated Gaussian b distribution with a lower 
A'h i-independent b c , Kim et al. (1997) concluded that b c 
over N Kl = 10 12 ' 



16 cm 2 increases as z decreases: 15 



km 



at z ~ 3.7, 17 km s 



at z ~ 3.3, 20 km s" 1 
2.3. This result 



at z ~ 2.9, and 22-24 km s" 1 at z 
has been explained by an additional heating due to the 
on-going He n rcionization, although the high b c value at 
z ~ 2.3 does not agree with any theoretical explanations 



(K im et a l. |1997| ; Haeh nelt fc Steinmetz |199J; Theuns et 
al. |2000b[ Schaye et al. |2000[ ). 

Subsequent studies on the z-evolution of b c have led 
to contradictory results. While b c is clearly dependent on 
AT H i (Kirkman & Tytler |1997| Zhang et al. |1997|), b c (N Hl ) 



and the mean b value at < z > = 2.7 (Kirkman & Tytler 



1997| ) and at < z > = 1.7 (Savaglio et al. |1999| ) does not 
show any noticeable difference compared with at z > 3. 
Combining the observations with the numerical simula- 
tions, Schaye et al. ( 2000| ) found that b c at the fixed over- 
density 5 — 0, bg—o, increases from z ~ 4.5 (6,5=0 ~ 14.5 



reionization at z ~ 3 and then decreases from z ~ 3 
(65=0 ~ 19.5 km s " 1 ) to z - 1.8 {b 5=0 ~ 14 km s" 1 ). 
Ricotti et al. ( 2000[ ) also found a similar increase in bs=o 
at z ~ 3, although their bs=o at z > 2.8 and at z < 2.8 
can be considered to be constant at bg = o ~ 14.5 km s _1 
and at b$ =0 ~ 21 km s _1 , respectively. On the other hand, 
adopting a slightly different approach for identifying ab- 
sorption lines instead of the Voigt profile fitting, McDon- 



ald et al. (2000) found that there is no b c evolution over 
2.1 < z < 4.4 at the slightly higher overdensity S = 0.4. 

Fig. ^l] shows the cutoff b values at the fixed column 
density Nm = 10 135 cm" 2 , 6 c ,i3.5j from the two power 
law fits in Sect. 4.3 as a function of z. Keep in mind 
that 6 c (A r Hi) from the iterative power law fit is an up- 
per limit, while 6 c (A r Hi) from the smoothed b distribution 
is a lower limit. The 6 c ,i3.5 value shows a slight increase 
with decreasing z from both power law fits, with a pos- 
sible local 6 Cj i3.5 maximum at z ~ 2.9 (with the caveat 
that the line list of Q0302-003 is generated by a different 
fitting program with respect to the other line lists. This 
could introduce an artificial result at z ~ 2.9 as shown 
in Sect. 5.4, although its redshift range suggests an influ- 
ence of additional heating if He 11 reionization does occur 
at z ~ 3). When all the values from both fits are av- 
eraged, 6 Cj i3.5 = 17.6 ± 1.6 at z ~ 3.75 is smaller than 
&c,i3.5 = 19.9 ± 1.2 at z ~ 2.1, but the difference is signif- 
icant only at the 1.44ct level. 

In simulations, b c is usually measured at a fixed over- 
density 6 rather than at a fixed column density. Translat- 
ing an overdensity into the corresponding column density 
is not trivial and depends on many uncertain parameters, 
such as the ionizing background and the reionization his- 
tory. If the simple law between S and A?hi by Dave et al. 
( 1999I ) 2 ! is assumed, then S becomes: 



A^hi - [0.05(1 + 5) 10° 



4z l 143 xl0 14 cm" 2 . 



(7) 



For <5 = 0, the corresponding A% 1 is 1.92 xlO 14 , 6.04xl0 13 , 
4.56 xlO 13 , 2.28 xlO 13 , 1.87 xlO 13 and 1.15 xlO 13 cm" 2 at 
<z> = 3.75, 2.87, 2.66, 2.13, 1.98 and 1.61, respectively. 
Although this conversion does not include the effects of the 
He 11 reionization, we assume that it is correct at least on 
a relative scale at z ~ 3.75 and at z ~ 2.1, where the Hen 
reionization would not affect the temperature of the IGM 
as strongly as at z ~ 3|^| Therefore, we only discuss the 



It should be noted that Dave et al. (1999) assume the QSO 



km s 1 ) to 



3 (6, 



'5=0 



19.5 km s x ) due to Hen 



dominated Haardt-Madau UV background without He 11 reion- 
ization. No other scaling laws between A'hi and 8 as a func- 
tion of z from simulations under different UV backgrounds are 
found in the literature. 

3 Note that Ricotti et al. ( |200C| ) found that 5 = corre- 
sponds to 7V H i = 5.62 x 10 12 cm" 2 at z = 2.85 with Hen 
reionization at z ~ 3, almost a factor of 10 smaller than the 
column density calculated by Eq. |^. This discrepancy indicates 
the importance of using the correct ionizing background in sim- 
ulations to constrain the temperature of the IGM. 
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relative behavior of b$=o as a function of z, in particular, 
our 6,5=0 at 2 ~ 2.1 with other results mentioned above. 

For the iterative power law fit (the power law fit to 
the smoothed b distribution), the b value at <5 = 0, b$=o , 
is 20.8 (17.2), 20.1 (17.2), 19.4 (17.8), 20.5 (17.5), 22.4 
(21.6) and 20.4 (18.8) km s" 1 at <z> = 3.75, 2.87, 2.66, 
2.13, 1.98 and 1.61, respectively. In the second panel of 
Fig. |l], bs=o is fairly constant with z as b$=o ~ 17-20 
km s _1 , with a possible local maximum bs=o ~ 22 km 
s" 1 at z ~ 2.9. The observed behavior of 6,5=0 is qualita- 
tively in agreement with the results from McDonald et al. 
( 200C ). While the observations agree with th e fairl y con- 
stant &,5 = q at z < 3 derived by Ricotti et al. (200C), they 
do not show the abrupt increase of 6,5=0 across z ~ 3 as 
large as ~ 7 km s _1 found by Ricotti et al. The observa- 
tions at z ~ 3.75 and at z ~ 2.1 agree with the results by 
Schaye et al. (200C) which show similar bs=o at z ~ 3.75 
and at z ~ 2.1. In addition, the observations do not show 
a strong decrease of 65=0 from z ~ 3 to z ~ 2 as large as 
~ 4 km s _1 as found by Schaye et al. (200C). However, 
note that, considering the large error bars of Schaye et 
al. ( 2000D , the significance of the decrease of 65=0 from 
z ~ 3 to z ~ 2 is not very strong and their result is 
not in disagreement with ours. It should also be recalled 
that we are using the scaling law between S and AThi es- 
timated from the QSO-dominated Haardt-Madau Thi- If 
this QSO-dominated UV background is underestimated at 
z < 2 as suggested by the evolution of the absorption line 
number density, the actual 6,5=0 at z < 2.4 can be higher 
than 6,5=0 in Fig. [|l]. 

The third panel of Fig. |T] shows the median b val- 
ues as a function of z measured for two column density 
ranges: N Hl = 1Q 13 - 1 - U cm- 2 and N Hl = 10 13 - 8 - 16 cm" 2 . 
It is rather difficult to interpret the z-dependence of the 
median b values. It could be constant at 1.5 < z < 4 with a 
small cosmic variance at z ~ 2. On the other hand, it could 
be decreasing with z at z < 3.1, if we discard the median 
b values at < z> — 1.61, based on a small number of ab- 
sorption lines. Although simulations correctly predict the 
shape of the observed b distributions, the predicted me- 
dian b values are typically 5-10 km s^ 1 smaller than the 



obser ved on es at all z (Br yan fc M achacek [20001 Machacek 
et al. [2000[ Theuns et al. EOOObp . 



The bottom panel of Fig. |2l] shows the power law slope 
of the A/hi - b c distribution, (Tt — 1), as a function of z 
(see Eq. |l|). No particular trend is apparent and (Tt — 1) 
shows very little evolution at z < 3.1. Note that the lower 
(Ft — 1) at < z > = 3.75 is in part due to the lack of 
lines with b < 15 km s" 1 and Ah, < 10 134 cm- 2 (Fig. |). 
When (Ft — 1) is averaged over all the measured values, 
(T T - 1) = 0.16 ± 0.03 at z < 3.1 is larger than (T T - 1) = 
0.07 ± 0.02 at <z> = 3.75. Due to the lower (T T - 1) at 
< z > — 3.75, it also seems clear from Fig. ^ that gas at 
lower overdensities is cooler at z < 3.1 than at higher z 
(keep in mind that a fixed column density corresponds to 
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Fig. 21. The z-evolution of 6 c ,i3.5, bs=o, median b and 
(Tt — 1) from the two power law fits in Sect. 4.3. Stars 
and triangles represent the parameters measured from the 
iterative power law fit and from the smoothed b distribu- 
tion, respectively. Vertical bars represent la errors. In the 
third panel, error bars in the x-axis represent the z ranges 
over which the median b values were estimated. Solid lines 
represent the median b values over A^hi = 10 131 ~ 14 cm~ 2 , 
while dotted lines represent the median b values over 
A/hi = 10 13 - 8 " 16 cm" 2 . 



a larger gas overdensity as z decreases due to the Hubble 
expansion) . 

If we assume Eq. g again and T - 59.2 b 2 for thermally 
broadened lines, (j T - 1) ~ 2.857 (Tt ~ 1)- When this 
simple conversion law is assumed, (jt — 1) = 0.46 ± 0.10 
at z < 3.1 and ( 7T - 1) = 0.18 ± 0.05 at z - 3.75. These 
(7t — 1) could be considered to be consistent with the 



results by McDonald et al. (2000) within their error bars, 
although their error bars at z = 3.9 ((7t— 1) = 0.42±0.45) 
and at z = 3 ((jt — 1) = 0.30 ± 0.30) are rather large. 
Their (jt — 1) = 0.5 ± 0.15 at z ~ 2.4 agrees with our 
(7T — 1) = 0.46. Our (77* — 1) is m argin ally in agreement 



with t he res ults from Ricotti et al. (200C) and from Schaye 
et al. ( 2000 ). Their (77- — 1) values are lower at z ~ 2 and 
higher at z ~ 3.7, but within the error bars. 



6.3. The optical depth distribution function 

While the one-point function of the flux is more closely 
related to observations, the one-point function of the op- 
tical depth is usually calculated in simulations (Zhang 
et al. 1998; Machacek et al. 2000). In the simulation by 
Machacek et al. (2000), the Hi optical depth r at which 
the maximum tP(t) occurs, r max , is r ~ 0.8 at z = 4, 
t ~ 0.09 at z = 3, and r ~ 0.013 at z = 2. Fig. M shows 



tP(t) as a function of r, which is calculated from the ob- 
served spectra at z < 2.4 and the spectra generated with 
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where all the symbols have their usual meanings in this 
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Fig. 22. The z-evolution of tP(t) as a function of r. The 
tP(t) values are calculated from the observed spectra at 
z < 2.4 and from the spectra generated with noise at 
z > 2.4. 

noise at z > 2.4. The r max value is r ~ 0.6 at < z > = 3.75, 
r ~ 0.13 at < z > = 2.87, and r ~ 0.08 at < z > = 2.66, 
respectively. The observed r max shows a behavior similar 
to the simulated results by Machacek et al. ( [2000 ). How- 
ever, r max converges to r ~ 0.04 at z < 2.4, not showing 
any z-dependence. Also there is no z-dependence of tP(t) 
at t < 0.04 and at z < 3.1. 

The optical depth r ~ 0.04 corresponds to F ~ 0.96, 
while r ~ 0.013 corresponds to F ~ 0.99, almost to the 
continuum level. As z decreases, the number of pixels with 
F = 0.96-0.99 increases. These pixels are noise-dominated 
by the limited S /N and the continuum fitting uncertainty. 
Therefore, instead of showing the expected z-dependence 
of r max and tP(t), the observed r max and tP(t) approach 
to an asymptotic r max ~ 0.04 value at z < 2.4 and an 
asymptotic rP(r) at r < 0.04 and at z < 3.1, respectively. 

At 0.1 < t < 3 (or 0.05 < F < 0.9), the observed 
tP(t) is simply a different way of viewing the one-point 
function of the flux. As z decreases, tP(t) decreases due 
to the expansion of the universe. At r > 4 (or F < 0.05), 
tP(t) starts to converge again since it typically samples 
saturated regions, again dominated by noise. 



study (Weinberg et al. |1997p . 

Table || lists the lower limits on at different z, 
together with the parameter values used for calculating 
f^ Hl . It should be noted that the temperature used to 
calculate Eq. || is smaller than the values derived in Sect. 
6.2, but in the present discussion we are only interested 
in the lower bounds on Of, . These lower limits on ^ Hl 
are consistent with flj, = 0.0125 hr 2 from the Big Bang 
nucleosynthesis analysis (Copi et al. 1995). These values 



also indicate that about 90% of all baryons reside in the 
Lya forest at 1.5 < z < 4. 

The lower bounds on the density parameter from the 

P ( F) 

one-point function, Q b , are given by Weinberg et al. 



( p97D as 

tf (F) > 0.021/1- 3 / 2 



' [^{~\nF) 1 IP-<P{F)dF/Ff-</ 2 ' 



1 + z 



H(z)/H a 
5.51 



0.70 

V10 4 K ) 



0.35 



1/2 



(9) 



where /3 7 = (2 — 0.77t) 1 = 1.6 — 1.8, where jt is the 
powe r law index of the equation of state (Weinberg et al. 
|1997D . 

Table || lists the lower bounds on along with the 

parameter values used to calculate il b . The D, b val- 
ues are larger than 0^ HI since thi is not the true mean 
H I opacity, but the effective opacity which underestimates 
the true opacity when absorption lines become saturated. 
The lower fif, limits from P(F) are about a factor of 
1.6 h - 5 larger than the Big Bang nucleosynthesis analy- 
sis, Q b = 0.0125 hr 2 . 

Our new lower bounds on fit, are a factor 
smaller than some of t he pre vious results, Q^, = 
0.03 /i~ 2 (Rauch et al. |l997t Zh ang et al. |l99S| ; 



of 1.5 
0.017- 
Burles 



et al. |199S| ; Kirkman et al. |2Q0C| ; McDonald et al. g000j), 
but still consistent with them within the error bars. How- 
ever, our lower fib bounds are not consistent with the 
derived n b = 0.005 /i" 2 -0.01 h~ 2 from the high D/Hmea- 



surements (Songaila et al. 1994; Rugers & Hogan 1996) 



6.4- The baryon density fif, 

We derived the baryon density, Qi,, from two proper- 
ties of the Lya forest, Tm (including all the Lya forest 
regardless of A?hi) and P(F) (only for the forest with 
Nm < 10 16 cm- 2 ). 

If thi < 1.24, the lower limits on the baryon density 
become 

nl Hl > 5.1 x 10 4 (1 - exp(-THi)) (i + z)- 2 - 5 (1 + n oZ ) Q - 2 ' 
xh-^Tg^r™, (8) 



7. Conclusions 

We have analyzed the properties of low column density 
Lya forest clouds (AT Hl = 10 12 5 " 16 cm- 2 ) toward 3 QSOs 
at 1.5 < z < 2.4, using high resolution (R — 45 000), 
high S/N (~ 25-40) VLT/UVES data. Combined with 
other high-resolution observations from the literature at 
z > 2.4, we have studied the evolution of the Lya forest 
at 1.5 < z < 4. Two parallel analyses have been applied 
to the datasets: the traditional Voigt profile fitting anal- 
ysis and a statistical measure of the transmitted flux. We 
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Table 9. The lower bounds on Qh 





THl 




o p ( p ). b 










(xh- 1 - 75 ) 


(xh- 15 ) 






1.61 


0.086 


> 0.010 


> 0.020 






1.98 


0.161 


> 0.016 


> 0.021 






2.13 


0.131 


> 0.012 


> 0.021 






2.66 


0.234 


> 0.014 


> 0.020 






2.87 


0.275 


> 0.015 


> 0.020 






3.75 


0.733 


> 0.013 


> 0.013 






a For 


T = 6 000 K, r = 


0.9 x 10~ 12 sec" 1 


(at 


z ~ 1.6), 


1.4 x KT 12 


sec (at z ' 


-2-3), 0.64 x 10" 12 


sec 


" x (at z ~ 



3.75) , and Cl = 1 (Haardt & Madau 1996; Hui & Gnedin 



1997 



For T = 
1.4 x 10 



Weinberg et al. |1997| ) 
6 000 K, I 



0.9 x 10" 12 sec" 1 (at 
(at z ~2-3), 0.64 x 10" 12 sec 



~ 1-6), 

(at z ~ 



3.75), Q =_L0 = 1-58, and H(z) = H n% (1+z)* ( Haard t 



& Madau 1996; Hui & Gnedin 1997; Weinberg et al. 1997) 



find that the general conclusions from both analyses are 
in good agreement. Although the results are limited by 
the relatively small number of lines of sight, we find the 
following properties and trends in the z-evolution of the 
Lya forest: 

1) The differential density distribution function of the 
lower column density forest (Am — 10 12 ' 5 " 14 cm" 2 ) does 
not evolve very strongly, /(iVm) iVg, , at 1.5 < z < 4, 
with (3 ~ 1.4-1.5 and with an indication of an increasing 
/? to j3 ~ 1.7 at z < 1.8. The higher column density forest 
(iVm > 10 14 cm -2 ) disappears rapidly with decreasing z. 
The observed slopes of f(Nm) for various column density 
ranges are considerably flatter than numerical predictions 
for the same column density thresholds. The same conclu- 
sions are drawn from the one-point function and two-point 
function of the flux. 

2) The line number density of the Lya forest with 
Am = 10 13 - 64 " 16 cm" 2 decreases continuously from z ~ 4 
to z ~ 1.5, dn/dz oc (1 + z ) 2 - 19±0 - 27 , without showing any 
flattening in dn/dz at 1.5 < z < 4. For the lower column 
density range at A/m = 10 131 " 14 cm" 2 , the number den- 
sity evolution becomes weaker: dn/dz oc (1 + z ) 110±0 - 21 . 
The line counts as a function of the filling factor show 
a negligible z-dependence. These results strongly suggest 
that the main drive in the evolution of the forest at z > 1.5 
is the Hubble expansion and that Th i changes more slowly 
than a QSO-dominated background at z < 2, suggesting a 
contribution from galaxies to the UV background at z < 2. 
When combined with the results from the HST QSO ab- 
sorption line key project at < z < 1.5 with the same 
column density threshold, there is evidence for a slope 
change in dN/dz at z ~ 1.2. 

3) Deriving the cutoff & c (A/hi) as a function of z de- 
pends strongly on the methods used and the number of 
available lines. However, the cutoff b parameter at the 



fixed column density Am = 10 13 5 cm -2 , 6 C ,13.5) shows 
a weak increase with decreasing z, with a possible local 
frc,i3.5 maximum at z ~ 2.9. Despite being substantially 
uncertain due to the uncertain conversion from the observ- 
able parameters b and Am to the theoretical parameters 
T and S, the cutoff b value at the mean gas density, bs=o, 
is fairly constant with z as bg—o ~ 18 km s , with a pos- 
sible local b c maximum at z ~ 2.9. The observed slopes 
of b c (Nm) do not show any well-defined z-dependence ex- 
cept for a flatter slope at < z > = 3.75 than at z < 3.1, 
possibly due to the severe blending of lines with low-6 and 
low- Am at higher z. 

4) The velocity correlation function and the step 
optical-depth correlation function confirm that stronger 
lines are more clustered than weaker lines and that the 
correlation strength increases as z decreases. 

5) The mean Hi opacity is well approximated by a 
single power law at 1.5 < z < 4, 7m oc (1 + z ) 3 - 43 ± 017 j 
without showing any flattening towards higher z. The sig- 
nificant scatter in the mean H I opacity is likely to be 
caused by continuum fitting uncertainty and by an inclu- 
sion/exclusion of high-column-density Hi systems (JVm > 
10 17 cm -2 ) towards different lines of sight. 

6) The lower limit on the baryon density derived both 
from the H I opacity and the one-point function of the flux, 
il^ ya ~ 0.013 /i" 175 , suggests that most baryons (over 
90%) reside in the forest at 1.5 < z < 4. The contribution 
to fib from the Lya forest does not change much with z 
at 1.5 < z < 4. 
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